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ABSTRACT 


Given the orbital elements of an earth satellite, the right ascension of 
its position, and certain S/C design characteristics, a number of other 
parameters can be computed which are useful in analyzing earth-S/C 
relations and in checking telemetered attitude parameters. For the 
case where an earth sensor is mounted at an angle to the S/C spin axis, 
the determination of whether or not the earth is in view by the sensor is 
of interest. 

For the program described in this paper, known orbital elements, a 
starting S/C position right ascension, S/C spin axis right ascension and 
declination, sensor-to- spin-axis mounting angle, angular width of earth 
sensor, and sun position are used to determine: 

• declination of S/ C 

• radius vector of S/ C 

• apparent angular radius of earth seen from S/C 

• angle from S/ C spin axis to S/C nadir (earth center) 

• angle subtended by S/ C sensor scan across earth 

• time to or from perigee 

• earth shadow relation to S/C position (when earth seen by S/C). 

The computations are made for a series of right ascensions at incre- 
ments of a chosen angle. Earth shadow calculations require information 
on sun position at time of S/ C position. This is updated from initial sun 
position data at time of S/ C perigee. 


iii 








CONTENTS 


Page 


ABSTRACT iii 

1. INTRODUCTION 1 

1.1 Aim 1 


1.2 Initial Parameters 


1.2.1 Orbit 1 

1.2.2 Attitude 4 

1.2.3 Spacecraft (S/C) 6 

1.2.4 Sun 11 

2. DERIVATIONS 13 

2.1 S/C Position • • • 13 

2.1.1 Right Ascension of S/C 13 

2.1.2 Declination of S/C 13 

2.1.4 Radius Vector of S/C 14 

2.2 Half Earth Angle 15 

2.3 Angle from S/ C Spin Axis to S/ C Nadir 15 

2.4 Earth Crossing Angle 17 

2.4.1 Scan Range 17 

2.4.2 Spherical Triangle 18 

2.4.3 Half Crossing Angle 19 

2.5 S/C Time To/ From Perigee 19 

2.5.1 Angular Approach 20 

2.5.2 Areal Approach 21 

2.6 Interpolation of Anti-Sun Line to Time Since Perigee 22 

2.7 Earth Shadow — S/C Relations 25 

2.7.1 Shadow/ Terminator 26 

2.7.2 Horizon/ Terminator 28 


v 



CONTENTS (Continued) 


Page 

3. PROGRAM 31 

3.1 APL 31 

3.1.1 Description 31 

3.1.2 Advantages 39 

3.1.3 Disadvantages 40 

3.2 Coding 41 

3.2.1 Program Design 41 

3.2.2 Initialization Routines 42 

3.2.3 Program Driver Routines 44 

3.2.4 Utility Routines 44 

3.2.5 Computation Routines 44 

4. APPLICATION 47 

4.1 IMP-6 S/C Parameters 47 

4.2 Orbital Elements 49 

4.3 Sun Parameters 49 

4.4 Nominal Attitude Values 49 

4.5 Comments on Results from IMP-6 49 

APPENDIX A — APL CHARACTERS AND CONVENTIONS 51 

APPENDIX B - APL CODED PROGRAM 55 

APPENDIX C - IMP-6: SELECTED RUN RESULTS 62 

REFERENCES 70 


vi 



FIGURES 


Figure Page 

1-1 Cartesian and Spherical Coordinates 2 

1-2 Orbit Referenced to Inertial Coordinates 3 

1-3 Orbit Ellipse Geometry 4 

1-4 Eulerian Angles 5 

1-5 Spin Axis/ Earth Sensor Mounting Angle 8 

1-6 Digital Solar Sensor Schematic 9 

1- 7 Cross Section of Earth Shadow 12 

2- 1 S/C Declination and Angle from Perigee. 14 

2-2 Half Earth Angle 15 

2-3 Angle from Spin Axis to Zenith Vector 16 

2-4 Scan Annulus Intersection of Earth 17 

2-5 Spherical Triangles for Determination of Earth Crossing 

Angle 19 

2-6 Differential Area of an Ellipse 21 

2-7 Perigee Right Ascension and Declination 23 

2-8 Update of Anti-Sun Line 24 

2-9 Anti-Sun Line and S/ C Angle 26 

2-10 Terminator Tangent and Anti-Sun Line 27 

2-11 Half Angular Diameter of Sun 27 

2-12 S/C in Shadow, at Terminator, and Sunlit 28 

2- 13 Geometry of Horizon or Terminator Detection 29 

3- 1 Block Diagram of Main Driver Routine 42 

3-2 Block Diagram of Limited Driver Routine 43 



SPACECRAFT ORBIT/EARTH SCAN DERIVATIONS, 
ASSOCIATED APL PROGRAM, 

AND APPLICATION TO IMP-6 


1. INTRODUCTION 

1.1 Aim 

The aim of this paper is to present a time shared, remote site, demand 
processed computer program, its derivation, and test results from it 
for analysis of selected orbit, attitude, and S/C parameter effects on 
earth sensor detections of the earth. The S/ C is assumed to be spin 
stabilized with an earth sensor mounted at an angle y to the spin axis 
and having a field of view or scan angle of S c degrees. 

For pre-launch analysis, the program may be used to simulate effects 
for nominal parameters which then would be useful in preparing attitude 
data processing programs. After launch, comparison of results from a 
simulation using estimated parameters and from computations on actual 
satellite data would produce deviations helpful in isolating problems. 

1.2 Initial Parameters 

The aspect — position and orientation — of a satellite can be uniquely 
determined by the specification of six parameters — three for position 
and three for orientation or attitude. This determination is essential 
for correct analysis of S/C telemetry data. In addition, certain S/C 
design characteristics affect the S/C sensor readings and subsequent 
telemetry readout. Finally, knowledge of the sun's position allows pre- 
diction of the S/C's view of the shadowed earth. 

1.2.1 Orbit 

S/C position is specified by reference to its orbit at a certain 
time. This reduces to three spatial coordinates — usually inertial 
geocentric Cartesian or spherical. Cartesian and spherical coor- 
dinates are related as shown in Figure 1-1. The Z axis is coin- 
cident with the earth's polar or spin axis, the X axis is in the 
equatorial plane intersecting the celestial sphere at the First 
Point of Aries, T , and the Y axis is defined to make a right ^ 
handed orthogonal coordinate system. An arbitrary vector R can 
be expressed in terms of the X, Y, Z components such that 
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Figure 1-1. Cartesian and Spherical Coordinates 


+ V + Zr 2 = R 2 - 

The vector R is also specified by the magnitude R, right ascen- 
sion A, and declination S where A is the counter-clockwise rota- 
tion of T to the R projection on the X-Y plane and S is the angle 
from there to the R vector. The transformation is: 

X = R cos S cos A 
Y = R cos 5 sin A 
Z = R sin 5 


The orbit is characterized by six independent orbital elements. 
The most common is the set of Keplerian elements consisting of: 

e — eccentricity or shape of orbit ellipse 
a — semi-major axis or size of orbit ellipse 
i — inclination of orbit plane to equatorial plane 
Cl — right ascension of ascending node 
0 p — argument of perigee or tilt of orbit 
t p — epoch or time of perigee. 
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Often included to define a particular point on the orbit at arbitrary 
time t is: 


2 7T 

M - mean anomaly - per ~^ x (t - t p ) 



Figure 1-2. Orbit Referenced to Inertial Coordinates 


Derived quantities such as semi-minor axis b = /a 2 - e 2 a 2 , 
perigee distance R p = a(l - e), apogee distance R a = a(l + e), 


and, from Kepler's third law, period T = 


4 7r J 
g 


a 3 (where g is 


earth's gravitational constant) are formed to simplify computa- 
tions and visualization. 


Because a satellite's orbit is not constant in time due to perturba- 
tions by the moon, sun, and non-spberieal nature of earth's gravi- 
tational field, the time rates of change of fi, 6 p , i, and T or fre- 
quent updating of elements are needed to provide accurate solutions 
of desired parameters. 
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Figure 1-3- Orbit Ellipse Geometry 


1.2.2 Attitude 

Attitude is the relationship of the S/ C coordinate system to the 
inertial (or other intermediate) coordinate system. This relation- 
ship is often expressed in terms of Eulerian angles such as roll, 
pitch, and azimuth (Figure 1-4): that is, a set of three rotations 
which transform one coordinate system of arbitrary orientation 
into a parallel orientation with the respective axes of a second 
coordinate system. 

Let: 


X, Y, Z = S/C coordinates 

X"', Y'", Z'" = inertial coordinates 

X', Y', Z' and X", Y", Z" = intermediate coordinates 
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Figure 1-4. Eulerian Angles 

roll angle 4> = rotation about X axis 

pitch angle v = rotation about new Y' axis 

azimuth angle 9 = rotation about final Z" axis 

The S/C coordinates can be expressed as a function of Euler angles 
and inertial coordinates in matrix notation by: 


X 


X'" 

Y 

= R ( 0 ) R (t?) R ( 6 ) 

Y ,w 

Z 


Z" 


where 



cos 6 sin 6 0 


cos rj 0 - sin rj 


i 

o 

0 

1 

R(0) = 

- sin 9 cos 6 0 

- R(t?) = 

0 10 

. R(0) = 

0 cos sin 0 


_ 0 0 1_ 


_sin rj 0 cos 7] _ 


_0 - sin cos 
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The inverse relation, solving for inertial coordinates, is: 


X"' 


V 

Y 1,1 

= R'(0) R'O?) R'(0) 

Y 

Z'" 


Z 


where the R' matrices are inverses of the respective R's. Since 
these matrices are orthogonal, the inverses are equal to their 
transposes e.g. 


R'(0) = 


cos 6 - sin 6 0 


sin 6 cos 6 0 


L o 


= [R (<9)] T = R (- 6 ) 


In many cases, the primary method is not by this approach to 
attitude solution but instead by the determination of the S/C spin 
axis in terms of inertial coordinates. Thence, evaluation of the 
S/ C spin rate and knowledge of the S/ C spin axis orientation 
effectively characterize the S/C attitude. This is equivalent to 
another set of Eulerian angles with respect to inertial and S/C 
coordinates: 

right ascension A. = rotation about Z'" axis 

co-declination (90 - 8) = rotation about new Y" axis. 

These two rotations transform Z'" axis in inertial coordinates to 
the Z axis of S/C coordinates. A final rotation derived from spin 
rate and a time since reference complete the attitude determination. 

1.2.3 Spacecraft 

For a spin stabilized S/C, the axis of rotation is coincident to the 
largest moment of inertia. These are often designated as the 
Z axis and I 3 , respectively. If the other principal moments of 
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inertia I x and I 2 are equal, then the S/C is dynamically balanced 
and the Z axis will tend to align with the S/ C momentum vector 
which is fixed in inertial space (when disturbing torques are 
absent). The effect of disturbance torques is slow such that ori- 
entation of the spin axis is approximately constant over the S/C 
period. A rapid reorientation may be commanded from the 
ground; in which case, the exact nature of the effect on the S/C 
must be known for accurate following of spin axis motion. 

To determine satisfactorily the Z direction, it is necessary to 
equip the S/ C with attitude sensors capable of measuring the 
direction of known objects from the S/C at selected times. These 
known objects may be the sun, the stars, the earth, and the 
earth's magnetic field. All of these are reasonably well defined 
as a function of time. For the S/ C discussed in this paper, two 
sensors are needed — an earth sensor and a digital solar sensor. 

The earth sensor is a narrow field of view instrument sensitive 
to the difference between the sunlit earth horizon and space and 
to the contrast between the earth terminator and the shadowed 
earth, the earth terminator being the boundary between the sunlit 
earth and the shadowed earth. When the S/C rotates, the sensor 
sweeps out an annulus on the celestial sphere with a radius the 
same as the sensor-to-spin-axis mounting angle. This angle y 
for particular spin orientation and orbit position is the major 
factor in whether the earth is scanned or not. Another factor is 
the field of view or scan angle S c . This is effectively the width 
of the annulus on the celestial sphere and determines the time 
the earth is seen. Pulses are initiated by photoelectric cell action 
whenever an horizon terminator is sensed thus providing time of 
earth and earth width measurements. 

For an earth horizon sensor for which the solution of spin axis 
equations have been based on detection of both horizons, it is 
important to determine if a terminator rather than an horizon is 
sensed. This is done analytically in most cases rather than 
through S/C electronics. 

The digital solar sensor provides a readout of the angle between 
the sun vector and the spin axis. With a wide field of view of 
about 180° centered at right angles to the spin axis,. it scans the 
entire celestial sphere every spin period, hence senses the sun 
each rotation (except when eclipsed by the earth). 
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SENSOR 
DIRECTION 
@ TIME t 


Figure 1-5. Spin Axis/Earth Sensor Mounting Angle 


The angle is determined by the response of a bank of photocells 
onto which a narrow slit allows sunlight to fall — the slit being 
perpendicular to the spin axis. The particular photocells excited 
represent a binary coding of the sun angle with the resolution of 
the sensor depending on the total number of unmasked cells. For 
instance, if the bank is a rectangular array of cells in rows and 
columns with the slit parallel to the rows (the columns parallel 
to the spin axis) and with the columns composed of different num- 
bers of alternating masked and unmasked cells, then the most 
significant digit is determined by the column with just two large 
cells — one of which is masked and cannot respond to the sunlight. 
In effect, this column determines which hemisphere the sun is in 
with respect to the spin axis. The next column has four cells each 
half the length of the ones in the first column. These two columns 
limit the sun angle value to a certain 45° sector. Each additional 
column further narrows the resolution on the sun angle until with 
eight columns the resolution is about 3/4°. This is due to the 
eighth column having 2 8 = 256 cells alternating 128 masked and 
128 unmasked such that 180° is divided into 256 bands as in 
Figure 1-6. 
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In addition to the narrow slit perpendicular to the spin axis, 
another narrow slit — the command slit — is parallel to the spin 
axis and supplies azimuth angle information. This slit has a 
180° field of view and excites at least one photocell in an addi- 
tional column of unmasked cells whenever the sun is sensed. 

The photocell signal is amplified and a centered sun pulse is 
generated from it. 

To complete the electronics necessary for interpreting pulse 
relations into attitude parameters, constant frequency counters 
are initiated or terminated by one of the various pulses supplied 
by the sensors. Thus a pulse from the digital solar sensor turns 
on a counter which is later turned off by a pulse from the earth 
sensor. The counter is read into the telemetry stream at the 
appropriate time and is converted on the ground into a relative 
time by: 


^•se/^se ^-se 

where C SE = counts between solar and earth sensor pulses, f SE = 
frequency of counter (counts/sec), and t SE = time between sun and 
earth sensing. 

Other relations are determined by counters which (1) count from 
sun pulse to sun pulse which provides spin period by: 

C /f = t 
'~sp / 1 sp l sp > 

(2) count from first earth sensor pulse from leading edge of earth's 
horizon to next earth sensor pulse initiated by trailing edge of 
earth's horizon which leads to scanned earth width by: 

(^ew/^ew) x (360°/ t sp ) = earth width angle, 

and (3) count from a certain bit of a certain word of the telemetry 
sequence to the sun pulse which gives a time relating the teleme- 
try to the attitude parameters by: 

C T s/fxS - *TS • 
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The sun sensor provides values for the spin axis which are 
located on an annulus around the sun direction with a radius 
equal to the sun angle. Similarly, the earth sensor horizon de- 
tections lead to the angle subtended by the S/ C spin axis and the 
S/C nadir or S/C to earth center vector. This angle in turn is 
the radius of another ann ulus centered at the earth to which the 
spin axis must point. The intersection of the two annuluses is 
the ambiguous solution of the spin axis when provided an earth 
sensor and a sun sensor. Previous knowledge of the attitude is 
usually sufficient to eliminate one solution. 

1.2.4 Sim 

The initial sun parameters are included in order to derive earth 
shadow and S/C relations. It is sufficient to know the right as- 
cension of the sun at the time of perigee, inclination (obliquity) 
of the ecliptic plane to the equatorial plane, and apparent angular 
diameter of the sun as seen from earth orbit. Specification of 
right ascension and inclination allow computation of declination 
and interpolation of sun position to time of S/C. Both computa- 
tions are discussed in Section 2. 

For a particular day and hour of the year, the American Ephemeris 
and Nautical Almanac supplies the right ascension of the sun, the 
inclination angle or mean obliquity of the ecliptic plane (which 
changes slowly and is approximately 23°27'), and the apparent 
angular diameter of the sun as seen from earth orbit (also slowly 
changing — near 32 minutes of arc). 

Since the sun has a finite distension and is not a point source, 
rays from different areas of the sun are tangent to the earth at 
slightly different angles. Those rays from the apparent edge of 
the sun directed tangent to the earth meet at a point producing a 
cone of shadow or umbra behind the earth. The points where the 
sun is only partially obscured by the earth is the penumbra. The 
umbra extends along the anti-sun line 1 348 640 km (Figure 1-7). 

When the S/C is in eclipse, its earth sensor cannot detect an 
earth horizon as it isn't sensitive to the differential between the 
shadowed earth and space. Location of the S/C with respect to 
the umbra is thus critical in knowing if attitude solutions from 
the sun and earth sensors will be possible at certain times. 
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Figure 1-7. Cross Section of Earth Shadow 


2. DERIVATIONS 


From the few initial parameters given in the previous section, a number of 
others can be derived which fully characterize the S/C's changing relation- 
ships to the earth and sun. Several relations are derived which could have 
been included as initial parameters; but where possible, derivations are 
made to simplify the input requirements. 

2.1 S/C Position 

In one orbit of the earth, a S/C revolves through 360° of right ascension. 
For each value in right ascension, a unique value of declination is found. 
For a series of right ascensions, a corresponding series of declinations 
is obtained, and for each S/C position pair of right ascension and declina- 
tion other parameters are derived. 

2.1.1 Right Ascension of S/C 

The series of right ascension inputs is determined by recognition 
of orbit sectors of interest. Thus, a preliminary general view of 
the orbit would require computations based on right ascensions 
10° apart over the entire 360° range. Brief analysis of results 
would indicate specific regions for further study. This could in- 
volve computations from RA's 1° apart over any 20° range, for 
example. To enable selected computation, specification of a 
starting RA, final or limiting RA, and RA increment is made 
when program is called. 

2.1.2 Declination of S/C 

For a particular S/C right ascension X., inclination of orbit j, and 
right ascension of ascending node Q, the S/C declination 8 and 
other angles of interest are found from spherical trigonometric 
laws of sines and cosines. Thus, from Figure 2-1, the interme- 
diate interior angle £ is calculated by: 


cos £ - - cos j cos 90° + sin j sin 90° cos (A. - fi) 

cos £ = sin j cos (\ - Q) . (1) 


The declination angle 8 is then given by: 
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Figure 2-1. S/C Declination and Angle from Perigee 


sin S - 


sin (X - fi) 

: z sin i 

sin C, J 


( 2 ) 


2.1.3 Radius Vector of S/C 

Also referring to Figure 2-1, the angular distance 9 0 from the 
ascending node to the S/C position follows from the derivation of 
8 by: 


cos - cos 8 cos (X - fi) + sin 8 sin (X - fi) cos 90° 

cos 0 Q = cos 8 cos (X - fi) . (3) 

The question of "which quadrant?" arises as & 0 can take on 
values from 0° to 360° while Eq. (3) limit results to 0°-180°. The 
desired results are obtained by setting: 


= 360° - 0 Oe 

when 

X 

- fi > 

180 

= Sn 

when 

X 

- n < 

180 


where 9 0g is derived from Eq. (3). 

The derivation of 6 0 is important in the determination of the 
angle from perigee 6 : 
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6 


( 4 ) 


where <9 p is the argument of perigee. This angle, in turn, is 
necessary to compute the radius vector to the S/C by an equiva- 
lent statement of Kepler's first law: 


:R P (1 + O 

1 + £ cos 6 


(5a) 


where R p is the distance from earth center to perigee and e is 
the orbit eccentricity. 

2.2 Half Earth Angle 

The apparent an gular radius p of the earth as seen from the S/C at It 
just depends on the earth radius R e and R. 

Re 

sin p - -jj- (6) 

Thus as the S/C orbits the earth, the angle subtended by the earth 
changes continuously, reaching a maximum at perigee — 2 sin -1 R e /R p 
— and a minimum at apogee — 2 sin -1 R e /R a . 



Figure 2-2. Half Earth Angle 

2.3 Angle from S/C Spin Axis to S/C Nadir 

The vector from the S/C to earth center is oppositely directed to the 
spacecraft zenith vector and as such is the same as the S/C nadir vector. 
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Figure 2-3. Angle from Spin Axis to Zenith Vector 


The angle from the spin axis to the S/C nadir can be determined if 
given the coordinates of the spin axis and the S/C nadir. This angle is 
necessary for comparison with the spin axis-to- earth sensor mounting 
angle y to indicate if the earth is viewed by the earth sensor. If the 
spin axis right ascension and declination are \g and Sg respectively 
and at time t the S/C position coordinates are \ and S z , then the angle 
a' between the two axes is found by: 


i jj - cos -1 cos (A.g - \ z ) (7) 

where the cosine functions serve to limit the value of \p to less than 
180° as required. Thence, 


cos a - cos (90° - Sg ) cos (90° - S z ) + sin (90° - Sg ) sin (90° - S z ) cos i p 
cos a! - sin Sg sin S z + cos Sg cos S z cos (8) 
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Since a' is the angle between the spin axis and the zenith vector, the 
angle from the spin axis to the nadir vector is: 


a = 180 ° - a' 


( 9 ) 


2.4 Earth Crossing Angle 

By the previous developments, the parameters are made available 
which are necessary for calculating the angle of the horizon-to-horizon 
crossing of the earth sensor scan. This angle is the same as the one 
derived from S/C telemetry from counts between the first earth horizon 
pulse and the second earth horizon pulse (assuming both horizons and 
not an horizon and a terminator are viewed). 


s/c 



Figure 2-4. Scan Annulus Intersection of Earth 


2.4.1 Scan Range 

Depending on whether a is less than y or a is greater than y, 
the intermediate angle v and subsequent angle comparisons differ. 

For a > y, define: 


v - r + y s c • 
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If 


T) < a - p (10a) 

then earth is not in scan range. 

For a < y, define: 


v - 7-^S c . 


If 

v > a + p , (10b) 

then earth is not in scan range. 

Where S c is the earth sensor scan angle, y is the mounting angle, 
and p and a are as developed above. If a = y, then the earth 
must be in scan range such that v = y in following development. 

.4.2 Spherical Triangle 

When the proper selection of intermediate steps leading to the 
crossing angle is made, the algorithm is simplified by holding 
true independently of the relation between a and y. 

Define a spherical triangle by the S/ C spin axis intersection of 
the celestial sphere (CS), the earth center vector (S/C nadir) in- 
tersection of the CS, the second horizon vector extended to CS, 
and the great circle arcs connecting these CS intersections. 

These three arcs are a, v, and p respectively as in Figures 2-5. 
By the law of cosines, the interior angle X is found from: 


cos p - cos v cos a + sin v sin a cos X 


or 


cos X 


cos p - cos v cos a 
sin v sin a 


( 11 ) 
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Figure 2-5. Spherical Triangles for Determination of Earth Crossing Angle 


2.4.3 Half Crossing Angle 

Finally, by the law of sines applied to the right spherical triangle 
containing 


sin g _ sin v 
sin X sin 90° 

or 

sin g - sin X sin v . (12) 

Twice g is the crossing angle from horizon to horizon of the 
earth sensor scan. 

2.5 S/C Time to/from Perigee 

Two methods, one angular and the other areal, can be used to solve for 
the time a S/C takes to move from perigee to any position in the orbit 
or to move from any position to perigee depending on which time is 
shorter. 
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2.5.1 Angular Approach 

This time is related to the area of a sector of the orbit ellipse 
by Kepler’s second law of planetary motion: the radius vector 
sweeps out equal areas in equal times. In addition, an angular 
relation to time is represented by Kepler's equation: 


M - E - e sin E 


(13a) 


where E is the eccentric anomaly, e is the orbit eccentricity, 
and the mean anomaly M is the product of the mean angular mo- 
tion or average orbital irate n = 27 t/T and the time t from perigee 
M = nt. Also, T is the period of the orbit. Thus, 


n 


(E - e sin E) 
t 


(13b) 


From the geometry of Figure 1-3 , E is related to the radius 
vector by: 


sin E 


R sin 9 
a ]/\ - e 2 


where 


= a (* ~ e2 ) 

1 + e cos 9 ’ 


such that 


sin E 


Vl ~ e 2 sin 9 
1+6 COS 9 


(14a) 


(5b) 


(14b) 


Substituting (14b) into (13b) and solving for t gives: 


t 


T (._■,( )/l - e 2 sin 6 \ / j/l - e 2 sin 6 \ 

2 t t l^ in \ 1 + e cos 9 ) e \ 1 + e cos 9 ) 


(15) 
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2.5.2 Areal Approach 

The form based on Kepler's second law is derived from the 
formula: 


— = T 
A 77 a b 


constant 


( 16 ) 


where t i§ the time to/from perigee, A is the area of the sector 
of the ellipse bounded by the perigee and position vectors, T is 
the orbit period and 77 a b is the total area of the ellipse. 

The sector area is determined by integrating from perigee to 
position the differential area. This also leads to the differential 
statement of Kepler's second law — the areal velocity is constant. 
This is true since the angular momentum of a S/C in orbit is a 
constant, 


h = R 2 6 - constant, 


and the differential area is 

dA = yR • Rd0 

as in Figure 2-6 which leads to the areal velocity 


(17) 


(18) 


d A 
d t 


1 „2 d& 

y R yy - constant. 


(17') 


The sector area is obtained on integrating (18): 



Figure 2-6. Differential Area of an Ellipse 
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de 


(19a) 



~ e i 

T R 2 d0 
Jo z 


a 2 ( 1 - e 2 ) 2 f e 

^ Jo ( 1 + e cos 0) 2 


where again 


a (1 ~ e 2 ) 

1 + e cos 9 


(5b) 


The integration is obtained from tables to give 


a 2 (l- e 2 ) 2 

A = ^ — 

A 2 


1 \ e sin 9 


(1- e 2 )l 1+ 6 cos5 (1 ~e 2 f 


tan' 


{]/r^ 


- e 9 

+ 6 tan T 


(19b) 


Then from (16), the solution for t is: 


t = — r A = 

7T3D 


T a 2 ( 1 ~ e 2 ) 
7rab 2 


e sin 9 
1 + e cos 9 


<l-e 2 f 


tan_1 (]/Fr 


- e 9 

- tsn- 




e sin 9 
1 + e cos 9 


— - tan- 2 (i/y^ 

( 1 - e 2 1 + 


- e 6 

7 tan 7 


( 20 ) 


Either Equation (15) or (20) for t relates the time to or from 
perigee of a S/C to the angle between the radius vector at perigee 
and the radius vector to the S/ C position. For purposes of this 
paper, (20) is the preferred form as it can be evaluated for 9 
angles up to 180° without requiring quadrant considerations. 

2.6 Interpolation of Anti-Sun Line to Time Since Perigee 

Due to the possibility of interpolating sun position information to the 
time desired from an initial sun position at a reference time, entry of 
a lengthy table of sun positions vs. time is avoided. The reference time 
is the time of the S/C at perigee just prior to the series to be computed 
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Figure 2-7. Perigee Right Ascension and Declination 

and the reference position is the right ascension coordinate of the anti- 
sun line at reference time, where the anti-sun line is oppositely directed 
to the sun position. For any series of computations, the initial action is 
to determine the right ascension of the perigee X. . This angle is a 
function of the argument of perigee 6 p , right ascension of ascending 
node n, inclination of orbit plane j, and perigee declination found from: 

sin 8 p = sin 9 p sin j . (21) 

The equatorial plane angle k from ascending node to perigee right 
ascension is: 

cos k = cos 9 p /cos S p . (22) 

Thence the R A of perigee is given by: 

\ = fl - * (9 p < fi) (23a) 

Again, quadrants must be considered to provide correct solution in all 
cases. Whence 

\ p = n + k (0 p > fi) (23b) 
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Figure 2-8. Update of Anti-Sun Line 


Additionally, the reference arc to the anti-sun line from the first point 
of Aries is determined by: 


cos REF - cos RAg cos S s (24a) 

with 

REF = 360° - REF of Eq. (24a) (24b) 

for RA S > 180°. 

After the S/ C declination through time computations have been made, 
the time r from reference perigee to specified orbit position is com- 
pared to the time t y the sun takes to apparently move through 360° of 
arc — the tropical year. This leads to the angle A through which the 
sun's radius vector moved when the S/C was at reference perigee to 
the time the S/ C was at the specified orbit position. This angle A is 
the increment to the reference anti- sun line position: 


A _ _r_ 
360° " t 


or 


A = (r x 360°)/t y 


(25) 
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where t y is the tropical year, 

t y = 365.24219879 - .00000614 x DT (26) 

and D T is the time since Jan. 1, 0 hr 0 min 0 sec , 1900 in centuries e.g. 
DT = .71 for Jan. 1, 1971. The angle A in the ecliptic plane is added 
to the reference arc KEF to define the new anti- sun arc. The com- 
ponents of this arc R Ag and 8 - of the new sun position are evaluated 
from the laws of sines and cosines: 

sin S s sin (REF + A) 
sin i sin 90° 


or 

sin 8g = sin i sin (REF + A) . (27) 

Thence, 

cos (REF + A) = cos RAg cos Sg + sin RAg sin Sg cos 90° 

or 

. cos (REF + A) 

cos RAg = ^ Qs § _ (when REF + A < 180°) (28a) 

Again, with quadrant considerations being made, 


cos RAg 


360 - cos (REF + A) 
cos Sg 


(when REF + A > 180°) 


(28b) 


2.7 Earth Shadow — S/C Relations 

To determine if the S/C is shadowed by the earth or if it senses a 
terminator rather than a shadowed horizon, it is necessary to compute 
the earth shadow or umbra and S/ C intersection. 
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2.7.1 Shadow/ Terminator 


Let S be the anti- sun direction and z be the unit vector of the S/C 
radius vector; then the arc /3 of the great circle connecting the 
two vectors is found from: 

cos fi = cos (90° - S z ) cos (90° - Sg) + sin (90° - S 2 ) sin (90° - S § ) cos /j. (29) 

where 

cos /j. - cos (RA z - RA§) . (30) 



Let vector t be the terminator tangent intersecting the orbit at 
P Q . Also, S £ is the angular diameter of the sun from earth orbit. 
Because of S E , tf will not be parallel to S but will be displaced 
toward S by 1/2 S E (Figure 2-11). 

Let the radius vector to the intersection of the orbit with the 

U A A 

terminator tangent be R Q and the arc between z (R 0 ) and S be /3 0 . 
The angle a 0 between R 0 and the earth radius R e at terminator 
tangent is then given by: 

Re 

cos cr 0 - -jT" - (sin p 0 of Eq. (6)) (31a) 
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Figure 2-10. Terminator Tangent and Anti-Sun Line 


The angle between S and R e is 90° - 1/2 Sj. since the angle be- 
tween T and "S is 1/2 S E and t is perpendicular to R e . If general 
angle cr is defined as: 

K 

cos a - (- sin p of Eq. (6)) (31b) 

for any orbit locus R and /3 is as defined in Eq. (29), then the 
conditions for S/C and shadow are: 

cr + P + \ Se 900 (32a) 



t 


Figure 2-11. Half Angular Diameter of Sun 
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Figure 2-12. S/C in Shadow, at Terminator, and Sunlit 


for S/C entering the earth shadow at terminator, 
a + + y Sg < 90° 

for S/C completely in shadow, 

CT + /3 + y Sg > 90° 


(32b) 


(32c) 


for S/C to scan at least one sunlit horizon. 

2.7.2 Horizon/ Terminator 

Finally, for the case of Eq. (32c), to determine if the S/C senses 
two horizons rather than an horizon and a terminator, compare 
the angle /3 between S and R with the angle a between R e and R 
where R e is earth radius to the terminator tangent, R is the radips 
vector to the S/C, and S is the anti- sun line direction. If a ter- 
minator and an horizon is sensed, 

P - a < 90° - y Sg (33a) 

and (32c) holds. 
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Figure 2-13. Geometry of Horizon or Terminator Detection 


If two horizons axe viewed, 

P ~ ^ > 90° - y Se 
and (32c) holds. 


(33b) 
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3.1 APL 

3.1.1 Description 

APL is a programming language used in conjunction with remote 
terminals connected by dataphones to a digital computer. The 
characters used for APL require a special "type ball" on the 
typewriter terminal. These characters are illustrated and de- 
fined in APPENDIX A. 

Many separate terminals can be connected to the computer at one 
time through time sharing. All individual work spaces are kept 
separate so that, in general, no interference from other users is 
felt. Use of the computer is through a demand processor rather 
than a batch processor of conventional card/ tape programs. This 
means access to the computer is in real time with no delay be- 
tween coding and response. 

The following is a discussion of only those aspects of APL utilized 
in the program of Section 3.2 and is not a comprehensive examina- 
tion of all of APL's facets. For that, see REFERENCES. 

Control of the computer from the terminal is afforded by a set of 
system commands initializing an individual's use of the computer, 
controlling and monitoring workspaces, and setting output and 
printout limits. Thus, after turning on the terminal, phoning com- 
puter, getting tone, cradling phone, and, depending on terminal, 
setting switches and awaiting lights, the "sign-on" command is 
given: 

'i nnnn : LLLLL 

where nnnn is an assigned, unique user number and LLLLL is a 
user-chosen lock word. The return key is struck on the 2741 ter- 
minal (return key, then ALTN key and 5 key struck simultaneously 
on the 1052 terminal) to let computer respond with a message. 
Note, all computer messages start on far left of paper, but the 
starting column for the user is indented six columns. 

After the computer message is printed, the workspace containing 
the program is loaded: 

) LOAD rf WWU 
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where WWWW is the workspace ID and the program can be run, 
modified, or erased. 

When finished with the desired modifications and the corrected 
workspace is preferred for future use over the old, the newly 
altered workspace can be saved: 

) 0,1 vn 

with the old workspace ID retained. If initial creation of the pro- 
gram had been performed after sign-on and no load command I 
made, then the new programming effort can be saved and assigned 
a workspace name (or new workspace ID assigned to current 
workspace) by: 

}qAT/7? rj TTyU 

where WWWW becomes the new workspace ID and can be loaded 
in future terminal sessions. 

To keep aware of the workspaces and their contents assigned to 
the user's number, several commands are available: 

) L TB 

lists names of workspaces in user's library. 

) 77 f!B 

lists user created function in the workspace currently loaded. 

'I '//[]?£ 

lists global variable names in that workspace, where global refers 
to all variable names except those used as arguments of functions 
or localized to function (as described later). 

When a different workspace is desired and it is necessary to re- 
move current contents without losing it for future use, use: 

) CLF.AH 

If any workspace is no longer necessary, it can be deleted by: 

T -7 \>J T-J *•? 

where WWWW is the workspace ID. If particular variable or 
function names of an active workspace are superfluous, use: 

) ERASE J'TJ KKKK LL rili*!*!*!-' ... 
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where the assorted letters represent a list of variable and func- 
tion names. 

When all activities at the terminal have been completed, sign-off 
is by: 

n t? t? 


where LLL is a word lock which must be entered with the user 
number when signing on. 

Several other system commands exist and should be familiar to 
anyone wishing to use the full scope of APL, but the preceding are 
sufficient for many purposes. See REFERENCES for sources with 
more detail and fuller coverage. 

The characters of APPENDIX A supply the basic tools for 
operating on values both singly (monadieally) and jointly (dyadically) . 
Using just these operators and numerical constants provides much 
more versatility and usefullness than any desk calculator. But the 
addition of the special symbols and their rules of association makes 
possible another order of activity. Definition of variables and 
programming of functions is thence possible. 

One APL convention is especially important to keep in mind. 

Whether in desk calculator mode or in programming mode, there 
is no hierarchy of character importance, e.g. x does not get 
priority over + . Instead, all computations are performed from 
the right of a line of symbols toward the left. Each operation 
made in turn. With judicious use of parenthesis, the exactly de- 
sired algorithm is created e.g. 

21-5+3*2x6 (1) 

is interpreted as 

21- ( 5+( 3* ( 2x6 ) ) ) (2) 

in APL, rather than as 

(21-5)+((3*2)x6) (3) 

in an hierarchial scheme. For results of (1) to be as for (3), the 
parenthesis must be supplied. 
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To associate a value with a variable name, the special symbol «- , 
the specification arrow, is used e.g. 

4«-lS0 

T-f-90 


or more involved, 


n«- a -L- 

At any time the value of a variable name may be checked by 
printing the name with no associated operations. 

c 

270 

where the 270 is the computer response. If a new variable is 
being created and its value is desired, the quad symbol H is used. 

135 

where 135 is the computer response and would be the response if 
E was requested. 

When a mistake of a character or characters on a line is made 
but before the return key is hit, correction of the error can be 
effected by backspacing to error and striking the ATTN key on 
the 2741 or the LINE FEED key on the 1052. The carriage skips 
a line (prints a caret - on the 2741) and deletes the error charac- 
ter and all characters to the right. Correct characters are typed 
in to finish line as desired. 

Using just variables and constants with the operators of APPEN- 
DIX A, the terminal can be used conveniently in desk calculator 
mode. A more powerful set of operators and the flexibility of 
variables makes it more useful than a desk calculator. However, 
the full impact of the versatility of APL comes in programming 
functions . 

Programming mode is entered by typing a del V followed by one 
of six different function header forms: 
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No Explicit Result 

Explicit Result 
With Dummy Variable A 


ry 77 W TJ 

X? A rpij 

No Arguments 

r -jprprj y 

P A Tprpp y 

One Argument 

vj v pfpfl y 

v /\ +x FTN v 

Two Arguments 


The arguments and result are local to the function, which means 
the calling sequence arguments and explicit result global variable 
have names independent of the function header arguments and ex- 
plicit result dummy variable names. Additional local variables 
can be designated for use in the function (but not available for use 
after each execution of the function) by including in the function 
header each additional local variable separated by semicolons; e.g. 

HP Til X ; U ; 7 ; '' 

where V n TP X is one of the header forms and ; U ; V ; are extra 
local variables. 

Simple examples of calling forms are given below for the six 
header forms: 

•p mp 

! P'j’i] theta 
! PA 11 CF FTP ALPHA 

FOP CE^TEx PASS 

p-t-fx ( PTP A ) 

TEi!P*-PP.BS FTP VOL 

After the function header, the next lines may include coding of 
mathematical expressions, definitions, function control, or formats 
for output. If a mathematical expression is formed and variables 
used outside the function are needed, it is not necessary to bring 
them into the function through arguments or statements analogous 
to FORTRAN'S COMMON. If not local to another function, they 
are global and thus available to any function or expression as last 
defined or evaluated. In addition, the result of any expression in 
a function becomes global unless it is the same as the dummy 
variable of or localized in the function header. The benefit of 
localizing na m es is in being able to use the same names in different 
functions without interference. 

After the function header is typed and the return key hit, the com- 
puter types 
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m 


on the next line, skips three spaces, and awaits typed input by the 
user. The next lines are typed in similar fashion, building the 
desired algorithm. 

If branching or iterations are required, the branch arrow - is 
used at the first of a|line followed by the destination of the branch. 
This destination may be the number of a line, an APL expression 
resulting in a line number, a line label, the number zero which 
causes branches out of the function, a blank space which branches 
out of the function and all levels of functions being executed, a 
nullivector such as i o which causes sequence to drop to next in- 
struction, and various combinations depending on values of vari- 
ables. The following illustrate the branch: 

-> 5 

+MA RK 

-*■(/!/?= 2 ) fV.A RK 

where the 0 compression of MARK produces the null vector so 
sequence falls to next line; while for AD = 2, the 1 compression of 
MARK results in MARK so next sequence is to label MARK. 

The label is attached to the first of a line and separated from the 
expression by a colon, e.g. 

T 51 r‘APK: X.T^A + BrC 

The purpose of the label is in not having to change the line number 
each time a change in the program adds or deletes line(s) previous 
to the given line. 

The use of alphabetic statements for output edit and other purposes 
is permitted by the use of a pair of single quote marks separated 
by any number of APL characters, e.g. 

[6] 'SOLUTION OR KFPLFP''S FPU ATT ON: 1 

For each execution of the function through statement [ 6] , the 
computer prints out: 

SOLI) TICK 0 R ¥ FPL F.P ' 0 EOUt 7. TON : 

If an evaluated expression is wanted on the same line, add a semi- 
colon after the last quote mark and include the expression: 

r 1 0 ~\ 'Ll ST A P CF TP A VF.L ED: * ;£>«- VFL x 7 
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This results in the expression being evaluated, associated with 
variable D, and printed without requiring a quad symbol: 

D IS TA N OF TP A VFL F.D: 3 2 0 

where VEL = 32 and T = 10. 

The appearance of double quotes alone on a line will cause a line 
to be skipped at execution time. 

The last character of the function must be the del — either at the 
end of an expression or by itself on the last line. This is done in 
order to leave programming mode and to return to desk calculator 
mode where the function can be tested. 

If errors in the function are detected or changes in it are desired, 
a very convenient set of function editing procedures is available. 
To produce an entire listing of the function and leave the terminal 
in programming mode ready for entry of another line, type: 

del, function name (without arguments), open bracket, quad, 
close bracket, e.g. 

vrarm 

This produces a function listing and next line number then awaits 
a new line to be typed. If addition of new lines is all that is re- 
quired, type them and strike a final del. If a previously typed line 
is seen to need correction, then type: 

open bracket, line number to be corrected, quad, approximate 
correction column, close bracket, e.g. 

r 19 ] ri2H9] 

computer types user types 

When return key is struck after last bracket, 

T 1 2 ] 32+XT*(AiC) 

@ 

is printed by computer, line is skipped, and ball is aligned at 
position designated by @. If BI was the variable name desired, 
the user back spaces one (or any number of spaces forward or 
backward necessary to align ball beneath mistake), types a slash / 
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to delete the two, types a one (or any number 1-9 or A, B, C, . . . 
for 5, 10, 15, . . .) to cause that number of spaces to be inserted 
in front of the - , and strikes the return key. The computer re- 
sponds by retyping the line, deleting all characters with a slash 
beneath, inserting spaces in front of all characters with a num- 
ber or letter beneath, and then backspaces to the first empty 
space to be filled. 

r 1 2 "! B 4-XI* ( Ai-C) 

@ 

User now types I and RETURN. The computer types next line 
number and awaits user action. 

If nothing more needs to be corrected, V is typed leaving the rest 
of the function intact. However, if a line needs to be deleted, type: 

[line number] e.g. 

r 13] m 

C OMPUTE R USER 

[ 6 ] 

On 2741, hit ATTN key. On 1052, hit LINE FEED key. Computer 
appropriately skips to next line, erasing the line from the function. 

r 7 ] 

If a line needs to be seen before deciding if it is to be deleted, 
type: 


[line number □ ] e.g. 

r 7 ] ran] 

r 8 1 PICK*- A FT* ( up < 1 o o ) 

rob @ 

If a new line is not desired, either go to another line or type V to 
leave programming mode. Otherwise type new line. 

To insert a line in place of another without printing old line, type: 

[line number] , expression e.g. 

f 8 ] r 10 J VALUE+PICKtPRTCF 

To insert a new line or lines between original lines, type: 
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[decimal number between two successive line numbers] e.g. 

mi r 1 3 . 5 ] 

ri3. 5] » NEE LINE DEE' 

[13.6] 'NEE LINE T”0' 

New decimal numbers are supplied in this case up to [ 1 3 . 9 ] after 
which r 1 4 ] is typed by computer. To provide more new lines for 
insertion, a smaller decimal line number should be supplied such 
as: 


r i3 . ii] 

The sequence new line numbers are formed is thence: 

r 1 3 . 1 1 ] , r 1 3 . i 2 ] , r 1 3 . 1 3 j , r 1 3 . 1 4 1 , .r 1 3 . 1 5 1 , r i ? . i o 1 

ri3.17],ri3.18],ri3.19],ri3.2], r 13.3],[13.4], 

[13.5],... 

After all corrections, deletions, and insertions have been made 
and the final V typed to leave programming mode, the computer 
renumbers all lines as integers; however, before the V is typed 
all line numbers are as they were last printed. 

These rules can best be understood by using them in various 
situations and with an aim to minimum effort and efficient repro- 
gramming. 

3.1.2 Advantages 

The choice of APL as the coding language for the derivations of 
this paper is based on several factors: 

• the APL characters and coding rules lend themselves to 
mathematical applications. 

• the development of the program is facilitated by use of a 
computer from a terminal in the manner of a desk 
calculator. 

• the capability of writing and preserving functions for 
performing various algorithms is available and simple 
to use. 

• the modification of preserved functions and variables is 
simple and complete. 
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• extensive knowledge of all of the language or of the 
computer actions is not necessary. 

• the making, detection, and correction of mistakes is 
immediate, and in general, nondistruetive of further 
activity at the terminal. 

• initializing functions for repeated loops is simple. 

• output from running a program is immediately received, 
checked, and rerun if necessary. 

• an algorithm can be debugged before being coded in 
another language. 

• arrays are handled easily in many complex applications 
by the use of the powerful operation set. 

3.1.3 Disadvantages 

Due to the relatively slow print rate of the terminal typewriters, 
a user encounters a lengthy delay in getting on a terminal when- 
ever the previous user has extensive printout. 

Whenever the computer goes down during program testing, much 
tedious retyping or program checkouts are required to get back 
to the pre-failure status. This is an especially frustrating problem 
when corrections have not yet been saved. 

All input parameters must be entered by hand on most terminals 
and/or tape input facilities are not easily utilized and are not gen- 
erally available from the terminals. If reruns are contemplated, 
semi-permanent storage of the parameters must be prepared with 
resultant loss of available storage. 

Formatting of output with column editing is to some extent incon- 
venient in that columns of figures must be carefully arranged 
rather than occurring as a natural result of the format designators. 
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3.2 Coding 


3.2.1 Program Design 

To correctly write programs is the obvious aim of any user of 
the APL programming mode. In addition, several other factors 
should be kept in mind while coding. The effect of these factors 
is to have an overall program design. 

If a number of similar conversions are made, it is more efficient 
and clearer to define a function which is initiated whenever the 
conversion is needed. A single function can often be rewritten as 
several functions in order to allow certain algorithms to be used 
independently of the others. In this case, control or driver func- 
tions are written which call various functions and supply inter- 
mediate and relating steps. 

The design of the spacecraft/ earth scan program resulted in 
several types of functions being developed to perform the various 
actions: 

parameter initialization routines, program drivers, utility 
(or trigonometric conversion) routines, and computation 
routines. 

For analysis of a desired orbit, the orbital elements, spin axis 
right ascension and declination, sun elements at time of S/ C 
perigee, and fraction of century since 1900 are related to the 
proper variable names in a parameter initialization function. 

This function can be named at later times allowing recomputation 
through the program without retyping the parameters. Several of 
these initialization routines are formed for the orbits of interest. 
There are two drivers: one references the entire set of deriva- 
tions of Section 2, the other is for S/C position through earth 
crossing angle computations only. The drivers set iterations 
through the routines for a series of S/C position right ascensions. 
The utility routines are designed to convert from APL character 
representation of trigonometric relations to a more recognizable 
and simple form. The computation routines are particular algo- 
rithms to compute one or several dependent parameters. 

The block diagram of Figure 3-1 illustrates the relationship of 
the functions to one another for the main driver routine TRUN. 
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REQUEST (CALL ) 
RESPONSE (RETURN) 


Figure 3-1. Block Diagram of Limited Driver Routine 


The simpler structure of the block diagram for the driver routine 
RUN is shown in Figure 3-2. 

3.2.2 Initialization Routines 

By simply typing an initialization function name, a set of orbit, 
S/C, and sun parameters are fed into appropriate global variable 
names which can thence be referenced by any other function. Re- 
placement by another set is a matter of typing another initializa- 
tion function name. The general form of the function is as follows: 

VP# where P# is the program name and 
# is the number of the orbit order 
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Figure 3-2. Block Diagram of Limited Driver Routine 


J - orbit inclination 

EPSLON - orbit eccentricity 

THE TAP - argument of perigee 

OMEGA *- right ascension of ascending node 

PER - period of orbit 

RA *“ radius vector at apogee 

RP - radius vector at perigee 

AS - (RA + RP) x .5 (computes semi-major axis) 

BS *“ (RA x RP) * .5 (computes semi-minor axis) 

RAL - right ascension of S/C spin axis 
DAL - declination of S/ C spin axis 

UVW (call to function to compute eccentricity relations) 
RAS - right ascension of anti- sun line at time of perigee 
DAS - ARC RAS (computes anti-sun line declination) 

DT 4 - fraction of century since 1900.0 at perigee 
SR semi-diameter of sun for orbit 
I *- inclination of ecliptic for orbit 

A routine PARAM to produce a partial listing of the parameters 
of P# is included to provide a check of the current parameters 
available for computation. 
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3.2.3 Program Driver Routines 

The drivers RUN and TRUN are used for repeated passes through 
function ORBIT for a series of input S/C right ascensions. Each 
successive RA is a constant increment of the preceding one. The 
increment, starting RA, and limiting RA are the arguments of the 
function. This requires three function arguments whereas the 
allowed function formats seem to allow but two. Actually, since 
the input arguments may be n dimensional vectors, each vector 
component may be a changing, function input. 

Thus for RUN and TRUN, the first argument on the left of the 
function name is the incremental value applied to successive 
RA's and the second argument is a vector with two components — 
the first is the beginning RA and the second is the limiting RA. 
For example, 

10 RUN 20 50 

produces three passes through function ORBIT for RA's of 20°, 
30°, and 40°. Note, 50° is not entered. 

The difference in RUN and TRUN is in the presence of the com- 
putation algorithms for time, sun, and shadow parameters in 
TRUN but not in RUN. Thus, RUN output parameters are limited 
to those from the function ORBIT. 

3.2.4 Utility Routines 

The utility routines for this program are all trigonometric. They 
were created to eliminate extra coding necessary to express in 
degrees the trigonometric relations of APL which are simplest 
for radian measure. The functional coding also resulted in more 
understandable reading of the functions. One radian measure 
trig function is also included for this reason. 

3.2.5 Computation Routines 

In initializing parameters, the function P# calls two computation 
functions. The first, UVW, produces three values dependent on 
the orbit eccentricity: 

U = 1 - e 2 
V = 1 + e 
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W = 1 - e 


for simplification of formulae used in other functions. The func- 
tion ARC computes for a given RA of the ecliptic plane the corre- 
sponding declination. 

Function ORBIT determines S/C declination and radius vector, 
half earth (radius) angle, S/C spin axis to earth center (S/C nadir) 
angle, and whether or not earth is in scan range and if so the 
crossing angle. A print statement of each parameter is also 
produced. 

Having explicit results are two functions for computation of time 
to or from perigee. TIME is the function for the basic equation: 


T 


area of sector 
total area 


x period. 


The area of sector depends on the angle Q from perigee and is 
computed from function AREA. The total area of orbit ellipse is 
a variable, set as 


7i x AS x BS. 

TRUN also performs computations as well as controlling program 
flow. It determines RA and S of S/C perigee, sets the sun refer- 
ence angle, and for each cycle through ORBIT computes RA and S 
of anti-sun line, angles between S, R, and R e , and terminator/ 
shadow bounds. 
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4. APPLICATION 

Use of the described program for computation of S/C orbit/ earth scan 
relations has been made for several satellites, though chiefly designed and 
analyzed for IMP-6 (IMP- 1). A description of that spacecraft and its early 
orbit parameters follows 

4.1 IMP-6 S/C Parameters 

IMP-6 was engineered with a digital solar sensor and an earth horizon 
detector. Combined with the aspect system electronics, the sensors 
form an optical aspect system which provides: 

• the angle between the S/ C spin axis and the sun direction, 

• the angle between the S/C spin axis and the local vertical, 

• the fraction of S/ C spin period between sun sensing and earth 
sensing. 

The IMP-6 S/C is a sixteen-sided drum, 72 inches in length, 53 inches 
across, with antennae and booms extending from the main body. It 
spins about the drum's axis of symmetry with its attitude sensors 
aligned on the same side of the S/C such that there is no angular bias 
to measurements relating the two sensor references. 

The digital solar sensor is as described in Section 1.2.3 except it has 
a nine bit encoder for the sun angle allowing more resolution — to 
about a third of a degree. 

A tenth bit supplies azimuth data from the command slit. New digital 
solar sensor information is read into the telemetry stream along with 
ten additional remainder bits every 81.92 seconds as the twenty bit OA 
scan. Repeat OA scan data is read out three times before new sun 
angle data is supplied. 

The earth horizon sensor is a small (1-1/2°) optical spectrum instru- 
ment mounted 90° from the spin axis. Thus, the mounting angle y for 
IMP-6 is 90° and the SCAN angle is 1-1/2°. Detection of horizons or 
horizon and terminator produces signals which are shaped into pulses. 

The aspect system electronics produce a sun centered pulse whenever 
the azimuth bit is set by detection of the sun in the command slit. This 
pulse initiates 1.6 kc counters which are on until 


47 



• the next sun pulse 


• the earth pulse produced by the next earth horizon sensing. 

These counters represent spin period and earth-in respectively. The 
spin period count relates to time by 

T OD = C CD (counts) x .000625 (secs/count) 

SP br 

where 1 / 1.6 kc = .000625 secs/count. The earth-in count gives the 
time from sun to first earth horizon in a similar fashion. 

In addition, two successive earth horizon pulses trigger and stop a 
counter which thereby supplies an earth-width count to the telemetry. 
This earth-width count provides the time the sensor takes to scan from 
an earth horizon to the next. When coupled with spin period, the scanned 
earth- width angle is produced. 

T 

EW 

— — x 360° = Earth-width angle (in degrees) 

*SP 

Finally, a pulse generated by the S/ C initiates a 1.6 kc counter which is 
then stopped by the next sun pulse. The resultant count is the sun-time 
count. The sun-time relates all of the sensor times to the S/ C teleme- 
try word times. 

The optical aspect counts are provided to the S/ C telemetry every 
81.92 seconds and 20.48 apart. The order in which the counts are pro- 
duced and are placed in the telemetry stream is 


• 

OAST - 

sun time count 

• 

OASP - 

spin period count 

• 

OAEI - 

earth- in count 

• 

OAEW - 

earth- width count. 


For IMP-6, the nominal spin period is about five seconds. 
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4.2 Orbital Elements 


The IMP-6 pre-launch prediction orbital elements were as follows: 

semi-major axis — 114 606.63 km 

eccentricity — .942 

inclination — 28.3° 

argument of perigee — 293.8° 

RA of ascending node — 171.9° 

The orbit actually achieved shortly after launch on 3/13/71 1616 hours 
had these elements: 

semi-major axis 
eccentricity 
inclination 
argument of perigee 
RA of ascending node 

4.3 Sun Parameters 

The sun parameters on 3/13/71 at 1700 hours were: 

apparent RA - 23 h 32 m 20.07 s 

apparent declination — -3° 32 , 49.5" 

semi-diameter — 16'06.98" 

obliquity - 23° 26 , 42L , 775. 

The anti- sun line is 180° or oppositely directed to the sun line. Thus, 
the declination is of opposite sign and the right ascension is 180° less. 

4.4 Nominal Attitude Values 

The nominal orientation of the IMP- 6 spin axis is perpendicular to the 
ecliptic plane. Specifically, the +Z axis is aligned with the south ecliptic 
pole. Thus, RA ~ 90° and S % -66.5° for spin axis. 

Looking "down" on the +Z axis, the S/C turns or spins in a counter- 
clockwise sense. The spin rate is about 11 rpm. 

4.5 Comments on Results from IMP-6 

For the first few weeks after IMP-6 was launched, two periods each 
orbit were favorable for scanning the earth. However, the scans near 


- 109 529.835 

- .9381564 

- 28.597879° 

- 301.97723° 

- 233.6023°. 
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orbit perigee were almost entirely in the earth's shadow. Thus, only 
one period each orbit provided earth scan data to the telemetry. The 
simulation rim of APPENDIX C reflects this situation as well as other 
range, duration, and angle characteristics of the earth scans. Improved 
correspondence between telemetry and simulation is possible by adjust- 
ment of the mounting angle, field of view, and spin axis parameters. 


50 



APPENDIX A — APL CHARACTERS AND CONVENTIONS 


APL KEYBOARD 
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Primitive Scalar Dyadic Functions a a b 


Definition or Example 


Primitive Scalar Monadic Functions e b 


Definition or Example 


*-* ”2*3.5 

1 

*~* 2*3.5 

■ 

*-* 2 + “3.5 

Add 

n 

No Change 

0*B *-► *B 


*--*■ +3.S 


4*1 .25 


1.75 ♦>* 3.5*2 


A raised to the power B *-► A*B 
2 *— 4*. 5 "3 -*-*■ "27*(1*3) 


(•fl)**A logarithm of B for base A 


Larger of A and B 

6.01 *~* 6 . 01 T 6.01 


-*■ A[B_ 

'3 *— "3r~7 


*-* 5 1 13 
— *■ 0 | 5 
«--*■ 14 | 5 
14 *-*■ 1 | 3.14 


Smaller of A and B *-*■ 4LB 
6.01 *--*■ 6.01L6.01 ”7 *--*■ "3 L - 7 


fl-(MJ*lfl*U ■»”* A | fl For A* 0 

B ■*-* A \B For 4 = 0. B20 
Domain Error *-*-4 10 For >1=0 , fl<o 
4 4 ■*■*• S | - 1 1 D.E. 0 |~4 


See Table of Mixed Functions 


See Table at right 


Relational Tests t 
Result is l if 
the relation is 
true) result is 
o if the relation 
is false. 


Subtract 


Multiply 


Divide 


Pcwer 


Logarithm 


Maximum 


Minimum 


Residue 


Change Sign 


Signum 


Reciprocal 


Exponential 


Natural 

Logarithm 


Magnitude 


A 

fl 

A*B 

Avfl 

A*B 

~UaB) 

A*B 

-Uvfi) 

0 

0 

0 

0 

1 

1 

1 

1 

0 

1 

0 

1 

1 

1 

0 

0 

1 

0 

0 

1 

1 

1 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 


Less 

Not Greater 
Equal 
Not Less 
Greater 
Not Equal 


0-B *-*• -B 


Sign of B: 1 


1*B *B 


(2.71828.. , J*B 

2.71828. . 

. -*— *• * 

4 *-*• *1.386294361. . . 

20.0855. . 

, *-*• * 


(2.71828. . . )S B 
1.386294361. . . *-*• S4 


S amH *~* •*(/ 
"0.693147. .. *■■*■ e.5 


Smallest Integer 2 fl *-*■ \B 
4 *--*• r 3 . 141 ”3 *— r"3.141 

101 — * 

noi 

Largest Integer i B «-«• LB 



3 *-*• L3.141 

4 *-*• l 3.141 

101 *-* 

1101 

Absolute value of B ** |B 



9.5*-*- |9.5 

9.5 *-* | 9.5 

o *-*■ |0 



Factorial fl*!B-i +-*■ IB For B21 , B an integer; i *-*■ J 
or Gamma Gamxna(B+lJ !B for non-integer B 6 *-*■ 
39916800 *--*■ ! 1 1 2.68344... *■-*• »2.3 

3.3283... *-*• ! “2 . 3 D.E. *-*• J _ 2 


Roll Random choice from ifl *--*■ ?b (Origin dependent) 

Random integer from 12345678*-*- ?s 


Pi times I B*3 .14159. . . *-*■ OB 6.283185... *•*• 02 


Not | Logical Negation 

| Table of Dyadic o Functions 


( l-fl*2 )* . 

.5 

*--*• 

0OB 

( 1-fl* 2 J* , 

. 5 


0O fl 

Arcs in 

B 

*-*■ 

(-l)OB 

Sine 

fl 


lOB 

Arc cos 

B 

*-* 

( -2 )OB 

Cosine 

S 

*-* 

2 OB 

Arc tan 

B 

*--*■ 

( - 3 ) OB 

Tangent 

fl 

*-*■ 

SOB 

( "l*B*2 ) * . 

, 5 

*--* 

( - 4 ) OB 

( 1+B*2>*. 

. 5 


40fl 

Arcsinh 

B 

*-•* 

(-5 JOB 

Sinh 

fl 


50B 

Arccosh 

B 

*-*• 

(-6 JOB 

Cosh 

B 


6oB 

Arctanh 

B 


(-7 JOB 
B in 

Tanh 

radians 

B 


7 OB 


6 space indent 


Keyboard unlocked 
at left margin 


Cn] line of function 
indent to some position 


Determining and Changing Terminal Input Hdde 


Input Mode 


Immediate Execution 


Evaluated Input 


□ Input 

Uu&atched quote mark 


Function Definition 
Line-edit Phase 1 
Line-edit Phase 2 


Any Input 


Any APL expression 


1 line of characters 
Any characters 


Line of function or 
editing oosaand 


Line modification 
by ‘visual fidelity" 


)L0AD weld 


» OB8UBST BS-backspace 
' (quote mark) 


V or ) CLEAR or )L0AD wsid 


V as rightmost character 
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Primitive Mixed Functions 


pp A *--* Rank of 4 
00S ♦♦ 2 


2 Spff ♦♦ 1 2 3 4 5 

67 0 9 10 


Dimension vector of 4 pS Entity vector pp4 ♦♦ Rank of 4 

pP ♦-* 4 pJ ♦♦ 3 4 p7 ♦♦ iO ppB ♦♦ 2 ppp4 ♦♦ 1 

Reshape 4 to dimension V B *— 3 4pil2 12p5 ♦♦ il2 2 Spf ♦-► 1 2 3 4 5 

7pP ♦♦ 2357235 OpP ♦♦ tO (lO)pP ♦♦2 6 7 0 9 10 


Make 4 into a vector ♦ ♦ (*/p4)p4 P ♦♦ ,P ,5 ♦♦ il2 p.7 ♦♦ 1 

Join two arrays along an existing coordinate. 4 and B conform if either p4 ♦♦ pS or (I*ipp4)/p4 ♦♦ pS 
or p4 *-* (I«ipp5)/pB or if 4 or 8 is scalar. A scalar argument is extended (reshaped) 
to conform to the other argument. 4111 

'f'.'JIS' ♦-* * TBIS ' ABCDABCD 4 2 2 2 111 

2 3 5 7,1 2 *-* 2 3 5 7 1 2 X ,X ♦♦ EFCBBFGB 4.MV.C1] 5 5 5 ♦♦ .4 3 3 3 where MM ♦♦ 2 2 2 

IJKLIJKL 4 5 5 5 333 

Join two arrays along a new coordinate. Two scalars laminate to a length- 2 vector. Two non-scalar 
arguments must have identical shapes. If one argument is scalar, it is extended to match the shape 
of the non-scalar argument. I is non-integral and shows what existing coordinate (s) the new coordinate 
lies next to. 4D 


'O' ,'K' «-♦ ’OX’ 
100.5 ♦♦ 100 5 


C.t.SJD ♦♦ 4BC 
DBF 

C, [1 .5 JO ♦♦ 

BE 

CF 

<7,(0. 5] 

'*’ *♦ 4SC where 

*•* 

C 'ABC' 

D ** 'DBF' 

PC 2] ♦♦ 3 


BLi 1] — 1 S 

9 




ABCD 

XClj] ♦♦123 

4 


X(1 3(3 2 1] 

♦♦ 3 

2 1 

'ABCDBFGBIJKL'LBl 

♦♦ BFCB 

Pt 4 3 2 1 1 3] ♦♦ 7 

5 

3 2 2 5 


11 

10 9 


IJKL 

First s integers 


iS ♦♦ 1 2 3 4 5 

i0 ♦♦ An 

empty 

vector 









5 12 5 


Pirst occurrence of 

4 

in V, or l+pV if no occurrence 



3 5 4 5 ♦♦ PiX 


P» 3 ♦-♦ 2 


5 6 5 » 5 ♦♦ 1 

PCPiS 2] ♦♦ 5 2 


5 5 5 5 



Take or drop the (|p[j]) first (if P£I]ao) or last (if P[J]<0) elements along coordinate I of 4 4BC 

"2 tP ** 5 7 2* * TBIS * «— • IS ' S*'TBIS' ♦♦ Empty vector _ 6tP ♦♦002357 2 3tX ♦♦ BFG 


The permutation which would order / (ascending — ii or descending — t) 

43.7 10 3.7 1.2 ♦♦ 4 l 3 2 t3.7 10 3,? 1.2 ♦♦2134 FC4K] gives ascending sort 


Compress 

<5 > 

Expand 

V/A 

V\A 

1010 /P~2 5 13 4 BCD 

10 1 0/X ♦♦S? 1 2 3 4 ♦♦ 1 0 1/[1]I ♦♦ 1 0 UB 10 11 1\X ♦♦ B FGB 

1 0 l\i 2 ♦♦102 9 11 9 10 11 12 / JKL 

Reverse 

<5 > 

Rotate 

♦4 

804 

Reflect 4 0i4 ♦♦ 4 3 2 1 DCBA IJKL BCDA 

4X ♦♦ BGFE ex — 0(1 }X ♦♦ EFGB 1 0 "lOX ♦♦ BFGB 

Circulate 4 left 8 places LKJI ABCD LIJK 


coordinate PCX] of result 
Transpose last two coordinates 


AB1 

2 1M — RX «-* BFJ 
COB 
DHL 


Membership 


Deal <3> 


Decode, 
Base value 


pVtl — pU 

XcP — 

P«i3 *♦ 1 1 0 0 



Mil ♦♦ Random selection of w elements from »y without replacement 
( t6)«276 ♦♦ vector with exactly two l’s 


52752 ♦♦ card deck shuffle 


Polynomial evaluation 


1011 7 7 6 ♦♦ 1776 
1614 0 10 ♦♦ 1034 


Radix conversion of 4 10 10 10 10T1776 ♦♦ 

to base B 16 16 16 16x1034 ♦♦ 


I 


XiC ♦♦ Value of polynomial with coefficients C 

(for several values of X if X is a matrix) 


76 0 1T3. 141592 ♦♦ 3 0.141592 (Split integer 

0 10 and fraction) 


Generalized inverse of 4 

4 has 

(pp4 )»2 and t/p4 

B minimizes Euclidean norm of error term 

8 has (ppB)tl 2 

4 

e 

5 


and 1+pB ♦♦ 1 t p4 105 97 

114 B 3 

9 

2 ♦♦ 2 9 5 

K3(pP)*.*0 l 2 ♦♦ Quadratic least-squares 


7 

10 

2 

fit of V 


<1» Restrictions on argument ranks are indicated byi S for scalar, V for vector, M for matrix, 4 or B for any. Except as 

the first argument of 8i4 or 5(4], a scalar may be used instead of a vector. A one-element array may replace any scalar. 

<2 > Arrays used 1 2 3 4 ABCD 

in exanplest P ♦♦ 2 3 5 7 X ♦♦ 5 6 7 8 X -* BFGB 

9 10 11 12 IJKL 

<3> Function result depends on index origin. 

<4> Elision of any index selects all along that coordinate. 

<5> The function is applied along the last coordinate; the symbols /, \ , and e are equivalent to /, \, and q, respectively, 
except that the function is applied along the first coordinate. If CS] appears after any of the sysfeols, the relevant 
coordinate is determined by the scalar 5. 

<6» Index origin affects left argument of % and index list of [4;...;4]. 


Composite Functions 


6/4 Reduction across last coordinate of 4. 

•/[S]4 Reduction across coordinate 5 of 4. 

4S.es Generalised Inner Product of 4 and B. 

A ».9B Generalized Outer Product of 4 and B. 

where S can be any of ♦-■»*er l | !o<s*» »*v 
B.G. +/V sums all elements of vector V. 

A * . x B is ordinary matrix product of 4 and 8. 
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APPENDIX B - APL CODED PROGRAM ROUTINES 


INITIALIZATION ROUTINES 
VP 1[P]V 

7 pi 

[1] J-*- 28. 597879 

[2] FPSLON+Q .9381564 

[3] Til F TA. P+- 301.97723 
[41 OMFCA*- 2 3 3.6023 
[51 PER+-S 012.525 

[61 RA-*-21 2 4 3 9.76 

[7] PP+-6 619 .9106 

T81 AS<-(RA+RP)*2 

[9] BF+iPAxRP)* 0.5 
[101 RAL+Q7.2 
[111 ZM7>“66.5 

[121 UVV 
[131 1^-23.445215 

[14] P.AS+-1 7 3.125 
[151 l?AF-*-APC RAF 

[16] DT+-0 .712 

[ 1 7 ] SP+- ( ( 16*60)+7*3600) 
v 

VP 2[ni'7 
T7 P 2 

[11 < T +- 2 8 . .925 

[2] F.PSIOF+O .939 560665 

[ 3 ] TUFT A P+ 3 0 2 . 5 3 3 

[4] OFF FA -<-2 3 3 . 19 6 

[5] PFR.-t-SQ 5 6.1 

[ 6 ] P.-1 "-210945.7 

[71 6619 .9106 

[81 AF+-(PA+PP)x 0.5 
[S] BS4-(PA*RP)*0 . 5 

[in] RA L+-8V . 9 

[111 ZM/>~64.8 

r 1 2 ] PF.V 

[13] J«- 2 3 . 4 4 5 1 7 3 3 3 

6 14] PAR*- 176.87 

[15] PAF+APC RAF 

[16 ] PT<-0 .7124 

[17] PP-f- (16*60)+5. 9*3600 
V 
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77*3 TP] 7 
7 . 7 ? 3 

r l] J*-2 8 . 8131 
f 2 1 EPS LOU*- 0.93745 
T 3 1 TJl^TA P+-3 02.7863 
[ 4 ] OilT'dA •‘-233.0821 

[5] PP.7+-5 H6 1.713 

[6] n ,4-<-2 1 1 0 1 2 . 6 

r 71 .-,.^6811.451 

rm AP+(PA+PP)* 0.5 

[9] ( .P/S x pp ) * 0 . 5 

no] /M/>8 9.3 
r 1 11 PAL<-~ 68. 5 
r 1 2 ] f/7 r - ; 

f 13 1 J*- 2 3.4551569 4 

r 14] ,<?•«- 1 7 8 

r 1 5 1 7.4 7^-/1 PC PA7 

f 1 6 ] PT*-0 .7125 

[174 7.7^-( I5ffi0) + 5 . 28^3 60 0 

C 

nn^rn-jn 
7 P 5 

[I] Jt-23 . 844 8 

6 2] EPS LOP*-') .93106 

[3] TPFTAPf-3 0 6.509 

[4] Oi’T <7/4-*- 2 29 .419 

[5] PFF+ 5973.352 
Tel P,4 -*-2 10 5 9 0.3 

r 7 ] i ? P «-.7 517.35 

[8] /!?<-(P/|+PP) xO . 5 

[9] BS*-( PA xRP ) * 0 . 5 

[10] PAL-*- 8 5 . 5 

[II] DAL+“ 69.8 
[12] UVP 

f 1 3 1 !*■ 2 3.44512222 

6 14] p/Cf-igo . 8 

[15] ZMP-MPP P.4 5 

[16] PT+Q .7126 

[17] 7 P*- ( 16v60) + l. 6 t3600 
7 


56 



3 r ";] n 

7 P 3 

[ 1] 7 <-2 8. 8181 

[ 2 ] FP.EL 9/;<-0 .93745 
8 3 7 THETA P<- 3 0 2.7283 
[4] OEECA+- 2 33.0 821 
8 5] PER*- 59 61.713 

7 6] RA+- 211012.6 
f 7] PP-*~Ci 811.451 

[8 ] AE-*-(.HA + PP) x0 . 5 

f 9 ] 9P<-(P/ : xPP)*0 . 5 

TIO] P/!7-^3 9 . 3 

8 11] P77<-~6 8 . 6 

T12] U VW 

f 131 J<-23. 45515694 

[14] P7l9<-17 8 

[15] DAS+ARC PAR 

[15] P7<-0 .7125 

r 17 1 S !?•*■( 15r60) + 5. 2 8 * 3 6 0 O 

7 

VP4[[]]? 

7 P4 

[1] 7^-29.9656 

[2] FPSLOT!-*- 9. 9 3 49 181 

83] 7.7 FTAP+-3 9 6.5338 

[4] OMEGA-*- 2 29.3 509 

[5] PER* 9585.09 
r6] ‘PA -<-21 12 86. 4 

[7] PP<- 710 6.718 

[8] 715-<-(P.4 + PP) xO . 5 

[9] 3S+(P.AxRP)* 0.5 

[10] P.4 7 -<-8 8 . 2 

[11] DA L<-~ 7 0.4 
ri2] PIP7 

[13] J <- 2 3 . 4 551825 

[14] P7.9-f-18 6.125 

[15] DAE <- ARC HAS 

816] 77-0.7125 

[17] 5P<- ( 16t60) + 3-:-3600 

n 
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7P5 r H] 7 
7 P5 

[I] J *-2 2 .8448 

[ 2 ] EPS LON *-2 .93106 
r 3] TJIETAP*- 3 0 6.509 

[4] OMEGA*-?. 29 .419 

[5] PEP*-59 7 3.362 

T6] PA*- 2 10 5 9 0 . 3 
[7] PP*-7 517. 35 

[S] AS+OfA + RP) x0 . 5 
r 9 ] RS*-(RA.*PP)*Q . 5 

[10] RA.L*-S 5 . R 

[II] P/1 69.8 

[12] UVW 

[13] I *-2 2 . 4 4512222 

[14] PA P«-19 0 . 0 

[15] DAS+APC PAP 

r 16] Pr-f-0.7126 

r 17] SF*-( 16*60 ) + l. 6*3600 

7 


VPAPAWniV 
7 PAP AM 

[1] ' OPPI T IP CLIP A TTOP : » 

[ 2 ] 'ECCENTRICITY: • '%FPSLON 

[3] 'PP./ISP. OF ASC. NOPE: '\o'<ECA 

[ 4 ] ' ARGUMENT OF PE PI CEE: ' 'TPFTA.P 

[ 5 ] ’ PFPICEF : * . pp 

[ 6 ? » APOGEE : ' \pa 

[ 7 ] » DECLINATION OF Z AXIS: ' ;DAL 

[8] 'FT. ASC. OF Z AXIS: ' \PAL 

r 9] » SEMIMAJOR AXIS: ' ;AS 

[ 10 ] 'SEMI MINOR AXIS: ' -BE 


[ 11 ] 1 R.A. OF ANTI-SUN LINE AT PE PTC EE TIME: ';PA, 

[ 12 ] . 'DEC. OF ANTI -SUN LINE AT PE PI CEE TIME: ' ; DA 

7 
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PROGRAM DRIVER ROUTINES 


SRUNlUlV 

V A? RUN N ; Y 

[i] y*/T[i] 

[ 2 ] START -.ORBIT 7 
[ 3 ] Y*Y+M 

[4] -+(.?<!![ 21) /START 

[5] 7 

V 

VTtfi/ffCD] v 

V A? TRUN tl-,Y-,LARP-,SP.D-,SIB 

[1] DEP*ASIN((SIN THETAP) xSIN J) 

[ 2 ] NR A P*A COS ((COS THETA P ) 4 COS DTP ) 

[ 3 ] RAP * ( ( 3 6 0 | OMEGA + NEAP) x ( THETAP< 1 8 0 ) ) + ( 3 6 0 I OMEGA - NPAP ) x ( THETA P > 1 8 0 ) 

[4] 7«-/?[l] 

[ 5 ] PEF*A COS ( ( COS PAS)» COS DAS ) 

[6] REF*(REF*RAS<1Q0 )+( 360 -REF) *FAS > 180 

[73 ' * 

[8] " 

[9] START : ORBIT 7 
no] Y*Y+M 

[11] RA Z«-3 60 | LAMB 

[12] LAR*((EAZ<RAP)x3G0)+RAZ-PAP 

[13] LARP*((LARf-lSa)xLAR) + (LAP>t<>0)x(3f>0-r J AP) 

[14] 1 TIME(DYS , HRS , MIPS) TO/PROM PERIGEE: • ; 365 24 60 T TIM* LAPP 

[15] *(N0T=1) /CON 

[16] TRP*TIME LAPP 

[17] TRP*( ( (P.AZ-RAP) SO ) x ( PKR-TPP) ) + ( {PAZ-RAP) >0 ) x(PEP+TRP) 

[18] TI!*TFP 

[19] TY*( 365 . 24219 879-6 . 1BE~6xDT) x24*60 

[20] PRS*(T!'iTY)x360 

[21] RDR*REF+DRS 

[22] DASB*ASI N ( (SIN PDP)xSIK I ) 

[23] P.ASB*A COS ( ( COS RDP)*COS MSS) 

[24] RASB*( (RDRS180)xPASB)*(1BO<PDP)x360-RASB 

[25] SPD*RASB, DASB 

[26] 1 RA AND DEC OF ANT I- SUN LINE: ' -,SRP 

[27] MU * A COS ( COS ( LAMB-RASB ) ) 

[28] UDAZ*90-DAZR 

[29] NDAS*9 O-DASB 

[30] BETA.*A.COS(((COS NDAZ)xCOS NDAS)+(SIN NDAZ)x(STp NDAS)xCOS MU) 

[31] SI C*A COS (PEAR) 

[32] SIB*B ETA ,SIO 

[33] *(90<SIG+BETA+SP) / OVTT 

[34] 'S/C IN SHADOW -- ANGLES *P.OM ANTI-SUN LIEF, R, AND RE: ' \STB 

[35] *con 

[36] 0t)TT:*(90>RFTA+SP-SIG) /TERM 

[37] 'ANGLES FROM ANTI-SUN LINE, P, AND RE: ' -,S.TB 

[38] *CON 

[39] TERM: ' TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: ’ ;STB 

[40] *COtl 

[41] CON: ' ' 

6 4 2] " 

[43] *(Y<tn 21) /START 

[44] " 

[45] " 

V 
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UTILITY ROUTINES 


vsinl mv 

V Y+-SIN A 
[1] Y^IocMtP 

V 

VCOoCQ] V 

V Y+-COS A 

[ 1 ] Y+200A+P 

V 

VT,4tf[n]V 

V X+TA N Y 
[1] X<r3OOYH80 

V 

W1 SJ/i/LQ] V 

V Y+ASIN A 

ri] y<-( "10,4 >xi80toi 

V 

VA COSllT]V 
v a 

Cl] y-*-( "20,4 )xl 80*01 

V 

v/ir/iprn] v 
v y^42vi// y 

Cl] Y*-( 3oy)xl80*ol 

V 

v4j?M//rmv 

V X^ ART AN Y 

Cl] y«- 3 oy 

V 
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COMPUTATION ROUTINES 


VARCUWV 

V 7*- ARC PA;Y 

[1] 7*-A COS ( ( COS RA)*SIN I) 

[2] X*-(ACOS(( COS I)*SIN Y ) ) x ( x ( l go -PA ) ) 

7 

wvw rmv 
v uvn 

Til U-*-1l -epslon* 2 

[2] V+\ +EPSLON 

[3] W+l-EPSLON 

7 

V0/?3I!T[]]V 

7 ORBIT RAZB ; CZ PDR RALPH -,CX-,CNU \SNU 
Til- NOT+O 

r 2 1 /lM)+36Ox(/MZS<0tfS(7/l) 

[3] 74Z7?-<-iMZ9 + 4fl0 

[4] LAMB*-PAZR 

r 5 1 LMOH+LA MB - OMEGA 

[ 6 ] C7+((SIN J)*COS ifWW) 

[7] Z+-ACOS C7, 

T 8 1 DAZB-ASIN(( (SIN J)*SIN I, MOM) iSIN Z) 
r g 1 THNOT-*-A.COS ( ( COS DAZB)*COS IMOM) 

[ 1 0 1 THHOT+-( TP. NOT* ( IMOM <180))+(360- THNOT ) x ( J.Wf»> 18 0) 

Till THETA*-THHOT- THETA.P 

T 1 2 1 R*-<RPx(.l + EPSLON))i(l + EPSLONxCOS TRET A ) 
f 1 3 7 PDP+RAZB , DA ZB ,F 

ri4l ' RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OR $/Cx ' \PDR 

r 1 5 ] RMO+ASIN REAP 

[16] PST->-ACOS(COS(. \RAL-RAZB)) 

[171 ALPHA-+t80-ACOR(({STt!(90-DAL))x<Srn(90-DAZn))x<.CO<HP$I))) + (COfn$0-l>AZ,))*(COS(90-nA7.R)'>'l 
[181 RALPH+RHO .ALPHA 

ri 9 l 'HALF EARTH ANCLE AND 5 UC Z AXIS TO FAP.TR CENTER ANCLE -. ' ; PALP!! 

7 201 -+(ALPHA>GA MM A ,ALPHA<GA MM A ) / ALPS TO AM , ALP IT GAM. 
r 21 7 'EARTH CROSSING ANGLE IS: * ; 2 xRHO 

r 2 2 ] " 

r 2 3 1 NV+GAMMA 
r 2 4 1 -+COMXI 

[25] ALPGTGAM: NU+-GAMMA + 0 , S*SCAN 
[261 ■* ( { ALPHA -RHO'i.NU) / OUT 

[27] +COHXI 

[28] ALPLTGAM: NU+-GAMMA- 0 . 5 xSCAN 

[29] -*•( (ALPHA+EHO ) <NU) / OUT 

[30] COMXI : CNUt-COS NU 
r 31 ] SNU+SIN NU 

[32] CX*-((COS RHO)-CNU*COS ALPHA) iSNUxSIN ALPHA 

[33] X+A COS CX 

[341 XI~ASIN(SNUxSIN X) 

[351 'EARTH CROSSING ANGLE-. ';2 *71 

[36] ” 

[37] *0 

[38] OUT-. ' EARTH NOT IN SCAN RANGE' 

[39] NOT+1 

[40] " 

7 

VTIMEimv 
7 T+TIME TH 
[11 TOTAR*-°ASxBS 
[2] T*-(ARFA TH)*PER*TOTAR 

V 

HARE A [ H] 7 
v A*-AREA THE j X 
r 1 ] X+U-EPSLONxCOS THE 

[2] 4+0. 5 x( 4 S* 2 )x(J/* 2 )x(( ( 2*U*1 . 5 ) x {ART AN ( ( ( Mi 10*0 . 5 ) xTAN( THE A 2 ) ) ) ) - ( ( EPS LON * U) * ( SIN THE) iX) ). 

7 
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APPENDIX C - IMP -6: SELECTED RUN RESULTS 


APPENDIX r -- IMP-6: SELECTED PUP RESULTS 


SELECT INITIAL PA RAtLETEPS : 
P 5 

PRINT SELECTED PARAMETERS : 


PAR AM 

ORBIT INCLINATION : 
ECCFNTBI CITY : 

PT.ASC. O v ASC. NODE: 
ARGUMENT CF PERIGEE: 
PERIGEE : 

APOGEE: 

DECLINATION OP Z AXIS: 
PT.ASC. OF Z AXIS: 
SEMIMAJOR AXIS: 

SEMI f UNO p AXIS: 

N.A. 0* ANfl-St'N LINE 

DEC. OF ANTI-SUN LINE 


AT 

AT 


29 . 8448 
0 .93106 
229 . 419 
306 .509 
7517.35 
210590 . 3 
”69 . 8 
85.6 

109053. 825 
39787.95033 
PEPTGFF TIME: 

PEPICEF TIME: 


.190.8 

" 4.64573413 


fl COMPUTE FOR s/C RIGHT ASCENSIONS 0 THROUGH 350 DEGREES 
n IN INCREMENTS 0* 10 DEGREES: 

10 TRUE 0 360 


RIGHT ASCENSION. DECLINATION , AND RADIUS VECTOR Qv $ /c . 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANCLE: 
EARTH NOT IN SCAN RANGE ' ' ’ 


0 23.54522847 210562.0841 

1.735819845 69.47493863 


TI ME (D V S , HRS , MINS ) TO/ FROM PERIGEE : 


29.17891525 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANCLE AND S/C Z AXIS TO EARTH CENTER ANGLE- 
EARTH NOT IN SCAN RANGE 


10 20.01769243 174535.3662 

2.094265077 76.0797296 


TIME (DYS, NRS , MINS ) TO/FROM PERIGEE: 


1 0 35.13592133 


RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOP OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH NOT IN SCAN RANGE 


20 15.73884165 113689.3357 

3.21607777 83.26526337 


TIME(DYS , HRS .MINS) TO/FROM PERIGEE: 


0 11 10.65838321 


RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


30 10.79971136 70790.66005 

5.169300334 90.90390251 

10 .33421005 


TIME (DYS, HRS, MINS ) TO / FROM PERIGEE : 

FA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P, AND PE: 


0 5 31.81984186 

194.3745392 "6.145003105 

163.8646718 84.83069967 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOP OF S/C: 40 5.364062118 46032.79237 

HALF EARTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 7.964353889 98.78925892 

EARTH NOT IN SCAN RANGE 


TIME(DYS , HRS , MINS) TO/FROM PERIGEE: 


0 3 5.513999278 
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RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR. OF S/C: SO 0.3333350251 31014.50401 

HALF EARTH ANCLE AND S/C Z AXIS TO EARTH CENTER ANCLF : 11.52824375 106.6316334 

EARTH NOT IN SCAN RANGE 

TIME(DYS ,1’P.S .HIES) TO/FROM PERIGEE: 0 1 55.32072163. 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR 0 F S/C: 60 "6.014154107 23564.20714 

HALF EARTH ANCLE AND S/C 7, AXIS TO EARTH CENTER ANGLE: 15.70424333 114.0803272 

EARTH NOT IN SCAN RANGE 

TIME (DYS.HRS, MIES ) TO/FROM PERIGEE: 0 1 17.76837369 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR 0” S/C: 70 “ll. 40205327 18300.4040 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 20.23258707 120.7678623 

EARTH NOT IN SCAN RANGE 

TINE ( DXS , HRS , ft I NS ) TO /FRO!' PERIGEE: 0 0 55.68119562 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 80 "16.27230278 15065.19630 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANCLE: 25.04758701 126.3506694 

EARTH NOT IN SCAN RANGE 

FINE ( P V S , HRS ,!'It!S ) TO /FRO}’. PERIGEE: 0 0 41.61265283 

RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 90 "20.4672577 12827.37668 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 29.81695059 130.5952763 

EARTH NOT IN SCAN RANGE 

TIl'E(DYS , HRS , MINS) TO/FROM PERIGEE: 0 9 32.02912672 


RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 100 "23.90438852 11272.51772 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANCLE: 34.45893871 133.3182512 

EARTH NOT IP SCAN RANGE 

TIRFfDYS , HRS , FIRS) TO/FRO!' PERIGEE: 0 0 25 . '11 509423 


EIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 110 "26.55407095 10159.46500 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 38.88837906 134.4938085 

EARTH NOT IN SCAN PANGS 

TIRE (D V S, HRS , FINS ) TO/FROU PERIGEE: 0 0 19.8716,2173 


RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR. OF S/C: 120 "28.41805652 9343.344212 

half EARTH ANCLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 43.05049365 134.2034511 

EARTH NOT IN SCAN RANGE 

TIFF (SYS, HRS , Mil’S ) TO /FRO!’ PERIGEE: 0 0 15.71535082 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR Of S/C: 130 ”29.5102508 8735.216665 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER. ANGLE: 46,9003774 1 32.6082967 

EARTH CROSSING ANGLE: 47.18668561 

TIFFi V V S , HRS , MINS) TO/FROM PERIGEE: 0 0 12.28609651 

PA AND DEC OF ANTI-SUN LINE: 194.5775538 "6.220259842 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 64.85954108 43.0996226 


RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 140 "29.84352844 8279.999004 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 50.38172354 129.8944486 

EARTH CROSSING ANGLE: 69.77085659 
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TIME ( DYS , HRS , MINS ) TO/FROM PERIGEE: 

RA AND DEC OF ANTI -SUN LINE : 

TERMINATOR SENSED -- ANGLES FROM ANTI -SUN LINE, R, AND 

00 9 .353978813 

194.5794176 “6.230032883 

RE: 56.37535809 39.61827646 

RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 

150 “29.42249642 7944.144221 

53.40565186 126.2269883 

86 .3047586 

TIME (DYS , HRS .MINS) TO /FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 

0 0 6.758035295 

194.5810672 6.230717078 

47.92961343 36.59434814 

RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANCLE: 
EARTH CROSSING ANGLE: 

160 “28.24138698 7708.7764 

55.83156054 121.727966 

98. 56258356 

TIME( DYS , HRS .MINS) TO /FROM PERIGEE: 

RA AND DEC OF 'ANTI-SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, P. , AND RE: 

0 0 4.385471208 

194.5825749 “6.231342394 

39.4319214 34.16843946 

RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE- 
EARTH CROSSING ANGLE: 

170 “26.28672807 7566.026542 

57.45834651 116.4767786 

107.0364784 

TIME ( DYS , HRS , MINS ) TO /FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 

0 0 2.149263481 

194.5839959 “6.231931768 

30.84341252 32.54165349 

RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 

180 “23.54522847 7517.355617 

58.04441161 110.5250614 

111.8345017 

TIME ( DYS , HRS , MINS ) TO /FROM PERIGEE : 

RA AND DEC OF ANTI-SUN LINE >> 

S’/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 

0 0 0.02250769569 

194.5853759 “6.232504155 

2 2. 28245 889 31.95558839 

RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 

190 “20.01769243 7573.164531 

57.37381169 103.9202704 

112.9391145 

.TIMFADYS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI -SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE., R, AND PE: 

0 0 2.19494566 

194.5867564 “6.233076715 

14.48642163 32.62618831 

RIGHT ASCENSION , DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 

200 “15.73884165 7753.212239 

55.35075994 96.73473663 

110 . 3452748 

TIME (DYS, HRS, MINS) TO /FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 

0 0 4.433223217 

194.5881788 “6.233666623 

10.8857525 34.64924006 

RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 

210 “10.79971136 80^7.438984 

52.05960223 89.09609749 

104 .117846 

TIME(DYS , HRS , MINS) TO/FROM PERIGEE : 

RA AND DEC OF ANTI-SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 

0 0 6.809482583 

.194.5896888 “6.234292892 

15.9043864 37.94039777 

» 





right ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE : 
EARTH CROSSING ANGLE: 


,-220 5.364062118 8516.897451 

47.7478373 81.21074108 

94.54896853 


TIME( DYS , HRS , MINS) TO /FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

S/C IN SHADOW -- ANGLES FROM ANTI-SUN LINE, P, AND RF: 


0 0 9.411168757 

194.5913421 6.234978569 

25.29119606 42.2521627 


RIGHT ASCENSION, DECLINATION , AND RADIUS VF.CTOP OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO FART’! CENTER ANGLE: 
EARTH CROSSING ANGLE: 


230 0.3333359251 9394.87012 

42.75764248 73.36836664 

80.63119201 


TINE (.DYS , HRS , MINS) TO/FRO M PERIGEE: 

RA AND DEC OF ANTI -SUN LINE: 

S/C IN SHADOW -- ANGLES FRO If ANTI-SUN LINE, P, AND RE: 


0 0 12.35269294 

194.5932113 ”6.235753804 

35.9439982 47.24235752 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF FARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: ■ 
EARTH CROSSING ANGLE: 


240 6.014154107 10489.05746 

37.4 506 5-866 6 5. 9196-7 278 

60 . 5081148 


TINE (DYS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 


0 0 15.79350703 

194.5953978 6.236660618 

46.94367487 52.54934134 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OR S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 

E APT H CROSSING ANGLE: 


250 11.40295327 11986.79861 

32.14739397 59.23213773 

24 . 2610309.8 


TIMZiDVS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P, AND PE: 


0 0 19.96818935 

1"94. 5980507 “6.237760828 

57.87078839 57.85260603 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH NOT IN SCAN RANGE 


260 16.27230278 14006.2818^ 

27.08931497 53.64033065 


TIME(DYS , HRS , MINS) TO/EROM PERIGEE: 


0 0 25.23948341 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 270 20.4672577 16717.8143 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 22.42780927 49.40472368 

EARTH NOT IN SCAN RANGE 


TIME(DYS , HRS , MINS) TO/FROM PERIGEE: 


32.19698918 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 280 23.90438852 20381.805 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH. CENTER ANGLE: 18.23615171 46.68174876 

EARTH NOT IN SCAN RANGE 


TIME (DYS, HRS , MINS ) TO / FROM PERIGEE : 


41. 85150603 


RIGHT ASCENSION, DECLINATION, AND RADIUS VECTOR OR S/C: 290 26.55407095 25416.63233 

HALE EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 14.53341198 45.50619154 

EARTH NOT IN SCAN RANGE 


TIME( DYS , HRS.,MINS ) TO/EROM PERIGEE: 


■ 0 0 56.04256214 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 300 28.41805652 32523.8843 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 11.30940521 45.79654889 

EARTH NOT IN SCAN RANGE 


TIME ( DYS , HRS , MINS ) TO /FROM PERIGEE: 


0 1 18.35600985 



RIGHT ASCENSION , DECLINATION , AND RADIOS VECTOR OF S/C: 310 29.5102508 42926.62579 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE-. 8.544815779 47.3917033 

EARTH NOT IN SCAN RANGE 

TIME (Dy S , HRS , MINS ) TO / FRORIPEPICEE : 0 1 56.36190792 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 320 29.84352844 58817.53312 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 6.225388085 50.10555138 

EARTH NOT IN SCAN RANGE 

TIME ( DYS , HRS .MINS ) TO/FROM PERIGEE: 0 3 7.555267506 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 330 29 . 42249642” ' 84063 . 20077 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 4.35141027 53.77301167 

EARTH NOT IN SCAN RANGE 

TIME{DYS, HRS, MINS) TO/FROM PERIGEE: 05 36.29487557 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 340 28.24138698 124186.2532 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 2.943987588 58.27203398 

EARTH NOT IN SCAN RANGE 

TIME(DYS , HRS .MINS) TO/FROM PERIGEE: 0 11 21.38627963 


RIGHT ASCENSION, DECLINATION, AND RADIUS VECTOR OF S/C: 350 26.28672807 178414.5863 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 2.048710405 63.52322143 

EARTH NOT IN SCAN RANGE 

TIME(DYS , HRS .MINS) TO/FROM. PERIGEE: 1 0 59.45282956 


SELECT NEW INITIAL PARAMETERS 
F2 

PRINT NEW PARAMETERS 
PARAM 

ORBIT INCLINATION: 28.825 

ECCENTRICITY: 0.939560665 

RT.ASC, OE ASC, NODE: 233.196 

ARGUMENT OE PERIGEE: 302.533 

PERIGEE: 6619.9106 

APOGEE: 210945.7 • 

DECLINATION OF Z >1X15: "64.8 

PT.ASC. OF Z AXIS: 87.9 

SEMIMAJOR AXIS: 108782.8053 

SEMIMINOR AXIS: 37368.9935 

R.A. OE_ ANTI-SUN LINE AT PERIGEE TIME: 176.87 

DEC. OF ANTJ-SUN LINE AT PERIGEE ' TIME : 1.356372094 


« COMPUTE FOR RA 20 THROUGH 40 DEGREES BY 1 DEGREE: 
1 TRUN 20 40 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OE S/C: 20 16.76778967 105959.973 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANC-LE: 3.450953903 83.81960339 

EARTH NOT IN SCAN RANGE 

TIME{DYS, HRS. MINS) TO/FROM PERIGEE: 0 9 22.71690996 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OE S/C: 21 16.34227441 100778.7751 

HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 3.628604704 84.58921916 

EARTR NOT IN SCAN RANGE 
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TIME(DYS ,HRS ,MINS) TO /FROM PERIGEE: 


0 8 39.86207636 


RIGHT ASCENSI Oil , DECLINATION , AND RADIUS VECTOR OF P/C: 
HAL? EARTH ANGLE AND. S/C 2 AXIS TO EARTH CENTER ANGLE: 
EARTH NOT IN SCAN RANGE 


22 15.91016496 95849.9314 

3.81546617 ^ 85.36306335 


TIME(DYS , MRS .MINS) TO/FROM PERIGEE: 


0.9080598329 


RIGHT ASCENSION,. DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


23 15.47156886 91171.10519 

4.01158548 86.14099344 

5.074560081 


TIME(DYS , HRS .MINS) TO /FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR. SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 


0 7 25.48327796 

180.3311163 " 0.143584952 

152.87018 85.98841452 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF-EARTH ANGLE AND 'S/C 2 ’AXIS TO ' EARTH CENTER ANGLE: 
EARTH CHOSSING ANGLE: 


24 15.02660031 86737.42785 

4 . 217094091 ” — 86\” 9928 5 - 925 * 
7.03730979 


TIME(DYS , HRS, MINS) TO/EROM PERIGEE: 0 6 53.24765124 

RA AND DEC OF ANTI-SUN LINE: 180.3513596 0.1523630937 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE , R, AND RE: 152.2990574 85.782996 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 25 14.57538033 82542.11343 

HALF EARTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 4.431756966 87.70850276 

EARTH CROSSING ANGLE: 8.312522715 


TIME{DYS, HRS, MINS) TO /FROM PERIGEE: 0 6 23 . 89147269 

RA AND DEC OF ANTI-SUN LINE: 180.3697947 " 0.1603571112 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P. , AND RE: 151.6982505 85 . 5682430 ? 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


26 14.11803678 78576.97534 

4.655873141 88.49775797 

9.190457175 


TIME (DYS , HRS , MINS ) TO/EROM PERIGEE: 0 5 57.13381319 

RA AND DEC OF ANTI-SUN LINE: 180.386598 " 0.1676435083 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P. AND RE: 151.069227 85.34412686 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


27 13.65470451 74832.85369 

4.889374522 89.29045065 

9.778368451 


TIMECDYS , HRS , MINS) TO/EROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED — ANGLES FROM ANTI-SUN LINE, R, AND PE: 


0 5 32.72064878 

180.401929 " 0.1742914617 

150.4134259 85.11062548 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EAHTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


28 13.18552537 ' 71299.96334 

5.132276031 90.08639827 

10.17784079 


TIMEiDYS, HRS .MINS) TO/EROM PERIGEE: 0 5 10.42284905 

RA AND DEC OP ANTI-SUN LINE: 180.4159316 " 0.1803633689 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 149.7322482 84.86772397 

RIGHT ASCENSION, DECLINATION, AND RADIUS VECTOR OF S/C: 29 12.71064927 67968.17304 

HALF EARTH ANGLE AND S/C 2 AXIS TO EARTH CENTER ANGLE: 5.384585236 90.88540985 

EARTH CROSSING ANGLE: 10.76593934 



TIME{ DXS , HRS , MINS ) TO /FROM PERIGEE: 0 4 50.03412732 

RA AND DEC OF ANTI-SUll LINE: 180.4287353 "0.1859154082 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P, AND RE: 149.0270501 84.61541476 


RIGHT ASCENSION , DECLINATION , AND RADIUS , VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


30 12.230229.17 64827.22585 

5.646302078 91.68728584 

11.1376163 


TIME(DYS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 


0 4 31.36902349 

180.4404567 "0.1909980783 

148.2991379 84.35369792 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: ' 


31 11.74443111 61866.90995 

5.917418628 92.49181B06 

11 . 31450627 


TIME(D V S , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P, AND RE: 


0 4 14.26096682 

180.4512002 "0.1956567457 

147.5497652 84.08258137 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR 0 * S/C: 
HALF EARTH ANCLE AND S/C 7, AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


32 11.25342417 59077.18862 

6.197918851 93.29878963 

11.30620679 


TIME(DXS, HRS, MINS) TO/FROM PERIGEE: 0 3 58.56044869 

RA AND DEC OF ANTI-SUN LINE: 180.4610599 ”0.1999321226 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 146.7801314 83.80208115 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


33 10.75738541 56448.29638 

6.487778409 94.10797491 

11.11293778 


TIME(DYS , HRS , MINS) TO/FROM PERIGEE : 0 3 44.13332221 

RA AND DEC OF ANTI-SUN LINE: 180.47011991 "0.2038607394 

TERMINATOR SENSED -- ANCLES FROM ANT, I -SUN. LINE , R , AND PE: 145.9913816 83.51222159 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH. ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


34 10.25649881 53970.80802 

6.786964472 94.91913953 

10.72555301 


TIME(DXS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI -SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P, AND RE: 


0 3 30.85923635 

180.4784558 "0.20747537 

145.1846069 83.21303553 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR 0 * S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


35 9.750955129 51635.68568 

7.095435567 95.73204036 

10.12290117 


TIMEiDXS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES From ANTI-SUN LINE, P, AND RE: 


0 3 18.63020563 

180.4861354 "0.2108054205 

144.3608461 82.90456443 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C 7, AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


36 9.24095185 49434.30864 

7.41314145 96.54642557 

9 . 264519515 


TIME (DYS, HRS, MINS) TO /FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P , AND RE: 


0 3 7.34931201 

180.4932196 "0.2138772834 

143.5210876 82.58685855 
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RIGHT ASCENSION, DECLINATION, AND RADIUS VECTOR OR S/C : 37 8.726692967 47358.48956 

H ALP EARTH ANCLE AND S/C Z AXIS TO EARTH CENTER ANCLE: 7.740023003 97.36203466 

EARTH CROSSING ANGLE: 8.071170407 

TIME {DYS , HRS , MINS ) TO/ FROM PERIGEE: 0 2 56.92953265 

RA AND DEC OF ANTI-SUN LINE: 180.4997631 “0.2167146561 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, P, AND PE: 142.6662713 82.259977 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR 0 f S/C: 
HALF EARTH ANGLE AND S/C 7, AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


38 8.208388835 45400.48023 

8.076012174 98.17859861 

6 .35980999 


TIME {DYS , HRS , MINS) TO/FROM PERIGEE: 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANCLES FROM ANTI-SUN 


0 2 47.29268563 

180.5058149 “0.2193388289 

LINE, R, AND RE: 141.7972914 81.92398783 


RIGHT ASCENSION, DECLINATION , AND RADIUS VECTOR OF S/C: 
HALF EARTH ANGLE AND S/C Z AXIS TO EARTH CENTER ANGLE: 
EARTH CROSSING ANGLE: 


39 7.686255962 43552.96934 

8.421031934 98.9958399 

3 .432678659 


TIME { DYS , HRS , MINS ) TO / FROM PERIGEE : 

RA AND DEC OF ANTI-SUN LINE: 

TERMINATOR SENSED -- ANGLES FROM ANTI-SUN LINE, R, AND RE: 


0 2 . 38.3684847 

180.5114191 0.2217689415 

140.9149989 81.57896807 
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